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Preface 
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This  thesis  Is  concerned,  with  the  measurement  of  ' Pu/  Pu  ra- 
tios in  an  isotopically  mixed  sample  by  using  high  resolution  alpha  and 
photon  spectroscopy.  It  demonstrates  an  attempt  to  use  alpha-gamma 
coincidence  spectroscopy  and  singles,  alpha  particle  and  L x-ray  spec- 
troscopy to  determine  these  ratios  for  low  activity  samples  of  less 
than  2000  DPM.  Coincidence  spectroscopy  proved  to  be  unfeasible  at 
these  activity  levels,  and  statistically  significant  results  were  not 
obtained  from  the  singles  method;  however,  procedures  were  developed  in 
this  experiment  that  should  allow  a continuation  of  the  experiment  to 
provide  improved  results. 

I wish  to  thank  several  people  for  their  contribution  to  this 
thesis.  My  sincere  gratitude  goes  to  Dr.  Richard  Hagee,  my  thesis  ad- 
visor, and  Dr.  George  John  with  whom  I discussed  my  progress,  my  dif- 
ficulties and  eventually  my  results,  and  to  Lt.  Col.  Savery  Stuckey  and 
Capt.  Ronald  Jefferies  of  McClellan  Central  Laboratory  for  their  in- 
terest in  this  project  and  support  In  supplying  isotopic  samples  on  a 
moments  notice.  I would  also  like  to  thank  Capt,  Piltingsrud  and  Capt. 
Stencel  of  the  USAF  Radiological  Health  Laboratory  in  verifying  a gamma 
spectrum  acquired  of  one  low-activity  sample  by  a run  in  their  low-back- 
ground chamber,  and  Mr.  Roland  Armani  of  Argonne  National  Laboratory  in 
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supplying  two  high-activity  Pu  samples  for  calibration  sources. 
Finally,  but  certainly  not  least,  a very  special  word  of  thanks  is  de- 
voted to  my  wife  Judy,  who  typed  this  thesis  and  whose  affection  and 
understanding  saw  me  through  this  project. 
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Abstract 

High-resolution  alpha  particle  and  photon  spectroscopy  is  applied 

239  240 

to  determining  the  ratios  of  Pu  and  Pu  in  isotopically  nixed  sam- 
ples of  less  than  a few  thousand  disintegrations  per  minute.  The  util- 
ity of  alpha-gamma  coincidence  spectroscopy  in  measuring  ratios  of 
239  240 

Pu  to  Pu  is  explored.  Experimental  equipment,  procedures  and  cal- 
ibration methods  are  discussed.  Minimum  activity  levels  necessary  for 
detection  and  quantification  are  presented  for  the  equipment  used  in  the 
experiment.  The  minimum  activity  level  for  quantification  for  an  eight- 

hour  acquisition  time  was  found  to  be  at  least  20,000  disintegrations  per 
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minute,  A method  of  determining  the  Pu/  Pu  ratio  is  also  pre- 
sented for  the  use  of  singles  alpha  particle  and  L x-ray  spectroscopy. 
This  method  is  based  upon  differences  in  uranium  characteristic  x-ray 
intensities  of  the  daughter  nuclides  of  these  two  isotopes  following  the 

internal  conversion  process.  A procedure  is  introduced  for  correcting 

241 

a complex  L x-ray  spectrum  for  the  presence  of  Am  in  the  sample  and 
its  contribution  of  neptunium  characteristic  L x rays  to  that  spectrum. 
Instrumental  and  statistical  uncertainties  prevented  reliable  deter- 
mination of  Isotopic  ratios.  Improvements  in  experimental  equipment  are 
recommended  to  reduce  uncertainties.  Major  radiation  Intensities  are 
tabulated  for  plutonium  isotopes  and  for  several  other  heavy  elements. 
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ISOTOPIC  ANALYSIS  OF  LOW-ACTIVITY  PLUTONIUM  SAMPLES 
USING  HIGH- RESOLUTION  ALPHA  AND  PHOTON  SPECTROSCOPY 

I.  Introduction 

In  the  fields  of  radiochemistry  and  health  physics  it  is  frequently 
desired  to  know  the  isotopic  composition  of  samples  containing  actinide 
series  elements.  Since  most  radioactive  actinides  decay  by  alpha  emis- 
sion, alpha  pulse-height  analysis  is  routinely  used  to  quantify  the 
amount  of  a given  nuclide  in  the  presence  of  others.  However,  there  are 

some  nuclides  that  emit  alpha  particles  with  energies  too  close  to  one 

239  240 

another  to  be  resolved  by  alpha  spectroscopy  alone.  Pu  and  Pu 
are  two  such  nuclides  that  are  often  found  together  in  varying  ratios 
and  for  which  there  is  no  simple  method  for  determining  that  ratio. 

Background 

239  240 

The  two  isotopes  of  Pu  and  Pu  emit  alpha  particles  with  en- 
ergies that  differ  by  only  a few  keV.  Pulse-height  analysis  techniques 
using  state  of  the  art  semiconductor  detectors  alone  are  inadequate  in 
resolving  these  two  isotopes,  but  instead  yield  only  a sum  count  of  the 
number  of  alpha  particles  emitted  by  both  isotopes.  The  ratio  can,  how- 
ever, be  measured  by  mass  spectrometry.  The  number  of  atoms  of  each 
isotope  present  may  be  determined  from  this  ratio  and  the  total  alpha 
particles  emitted,  (determined  by  a spectrometer  with  known  efficiency). 
While  this  procedure  has  proven  reliable,  it  requires  additional  time 
and  effort  to  prepare  separate  samples  for  both  the  alpha  counting  system 
and  the  mass  spectrometer.  If  only  one  sample  is  present,  recovery  of 
the  sample  after  mass  spectrometric  analysis  is  difficult.  The  overall 
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efficiency  of  the  procedure  for  detecting  and  quantifying  Pu  ratios 
could  be  improved  if  the  procedure  could  be  reduced  to  a single  meas- 
urement on  a single  sample. 

Problem  Statement 

The  purpose  of  this  thesis  was  to  explore  and  determine  the  feasi- 
bility of  using  high-resolution  alpha  and  photon  spectroscopy  together 

239 

to  resolve  the  isotopes  of  heavy  elements  in  general,  and  those  of  Pu 
and  240Pu  in  particular.  The  object  was  to  find  a method  of  determining 
these  Isotopic  ratios  by  a relatively  simple  method  without  the  use  of  a 
mass  spectrometer. 

Scope 

The  study  was  primarily  experimental  in  nature  and  was  restricted  to 
applications  of  currently-available  alpha  and  photon  detection  systems. 
The  study  also  primarily  addressed  samples  having  low  levels  of  activity. 
The  activity  levels  varied  from  a few-tens  to  several  thousand  disinte- 
grations per  minute  and  were  the  guide  for  determining  practical  sensi- 
tivity levels  for  detection. 

General  Approach 

Two  approaches  were  studied j first,  the  techniques  of  coincidence 
spectroscopy  were  examined  for  applicability  to  the  problem.  Second,  the 
information  contained  in  alpha  and  photon  spectra  were  examined  to  see  if 
the  information  contained  in  each,  when  taken  together,  might  provide 
sufficient  information  to  calculate  the  isotopic  composition  of  the  sam- 
ple 

Coincidence  spectroscopy  uses  the  output  pulses  of  a coincidence 
circuit  to  turn  on,  or  "gate",  a multichannel  analyzer  for  storage  of  a 
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pulse  only  when  that  pulse  Is  In  coincidence  with  a second  pulse  meeting 
certain  criteria.  Specifically,  an  alpha  spectrum  was  acquired  only 
when  gamma  radiation  emitted  by  one  of  the  Isotopes  present  was  simulta- 
neously detected  by  a second  detector,  * 

The  second  approach  was  based  upon  the  fact  that  although  both  Pu 
isotopes  emit  alpha  particles  and  characteristic  x rays (following  in- 
ternal conversion) of  the  same  energy,  the  difference  in  their  x-ray 
emission  rates  per  alpha  disintegration  is  sufficiently  different  that 
the  disintegration  rate  of  each  isotope  could  be  calculated. 

Although  the  first  approach  proved  to  be  negative  at  the  activity 

levels  of  interest,  the  second  approach  appears  to  offer  some  promise. 

Using  the  procedures  developed  in  this  work,  and  perhaps  a different 
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detection  system,  one  should  be  able  to  determine  the  Pu/  Pu  ratio 
in  a single  measurement. 

The  remainder  of  this  report  is  devoted  to  a detailed  description 
of  equipment  and  procedures  used  in  this  analysis.  Chapter  II  provides 
the  theory  for  the  two  approaches  pursued.  Chapter  III  gives  a complete 
description  of  the  experimental  equipment  and  electronics.  Chapter  IV 
contains  experimental  procedures  and  calibration  methods.  Chapter  V dis- 
cusses the  methods  of  data  reduction  and  error  analysis,  and  Chapter  VI 
presents  results,  conclusions,  and  recommendations  for  future  work. 
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( II.  Theory 


Radioactive  Decay  Process  1 

239  240 

As  Is  the  case  with  many  radioactive  actinides,  Pu  and  Pu 
decay  by  alpha  emission  to  either  ground  or  excited  states  of  the 
daughter  nuclei  JJU  and  J U respectively. 

The  simplified  decay  schemes  depicting  major  branching  of  these 
two  nuclides  are  shown  in  Figure  1.  It  must  be  noted  that  the  decay 
schemes  of  these  two  nuclides  are  much  more  complex  than  shown  with 
many  higher  excited  levels  of  the  uranium  being  infrequently  populated 
by  emission  of  lower  energy  alpha  particles.  (Ref  1:427,429). 


Figure  1.  Decay  Schemes  for  Plutonium  Isotopes 


Whenever  a nucleus  is  formed  in  an  excited  state  for  which  the  exci- 
tation energy  is  insufficient  for  nuclear  particle  emission,  the  mode  of 
deexcitation  is  by  electromagnetic  transition  or  internal  conversion. 

The  electromagnetic  transition  consists  of  the  emission  of  a gamma  ray. 
The  internal  conversion  process  is  that  of  ejecting  an  orbital  electron. 
The  two  processes  are  competitive  and  the  ratio  of  internal  conversion  to 
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gamma  emission  is  known  as  the  Internal  conversion  coefficient.  This  co- 
efficient can  be  determined  by  rigorus  quantum  mechanical  calculations, 
or  determined  by  experimental  measurement  of  internal  conversion  electron- 
to-gamma  emission  ratios.  The  conversion  coefficient  is  dependent  upon 
several  factors.  As  the  energy  of  the  excited  state  increases,  the  in- 
ternal conversion  coefficient  rapidly  decreases.  The  conversion  coeffi- 
cient is  strongly  dependent  upon  the  parity  and  spin  differences  of  the 

excited  and  the  ground  state.  Furthermore,  the  coefficient  Increases 

3 

rapidly  with  atomic  number,  approximately  as  Z over  much  of  the  range. 

In  general,  internal  conversion  is  found  most  often  in  heavy  elements 
when  low  energy  gamma  rays  are  emitted,  (Ref  2:428). 

A vacancy  produced  in  an  atomic  shell  from  internal  conversion  causes 
electron  transitions  from  outer  shells,  a cascade  effect,  where  the  law 
of  energy  conservation  may  be  either  satisfied  by  the  simultaneous  emis- 
sion of  a photon  (characteristic  x ray)  or  by  the  ejection  of  an  electron 
from  an  outer  shell  (Auger  electron).  The  x-ray  emission  and  Auger  tran- 
sitions are  also  competitive,  but  x-ray  emission  is  favored  for  elements 
of  high  atomic  number. 

The  process  of  internal  conversion  follows  the  decay  of  the  parent 
nuclide  to  the  excited  state  of  the  daughter  nuclide,  thus  in  the  case  of 
plutonium  alpha  decay,  x-ray  emission  is  characteristic  of  the  uranium 
nucleus.  Although  the  process  does  follow  alpha  decay,  the  emissions  of 
alpha  particle  and  x rays  are  essentially  simultaneous. 

The  lower  excited  energy  states  of  uranium  are  highly  internally  con- 
verted | giving  rise  to  many  more  conversion  electrons  and  characteristic 
x rays  than  emitted  gamma  rays.  A summary  of  major  radiation  abundances 
for  some  plutonium  Isotopes,  along  with  those  of  several  other  heavy 
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elements,  are  tabulated  In  Appendix A>  The  most  abundant  gamma  ray  of 
239 

J Pu  is  that  of  the  51.6  keV  transition  which  is  present  in  an  abun- 

-4.  240 

dance  of  2x10  gamma  rays  per  disintegration.  In  the  case  of  Pu, 

the  most  prominent  gamma  ray,  45  keV,  is  only  slightly  more  abundant 
with  4,5x10  gamma  rays  per  disintegration.  In  both  of  these  pluto- 
nium isotopes,  the  most  abundant  photons  are  the  L-series  characteristic 
x rays  from  internal  conversion  which  are  approximately  100  times  more 
abundant  than  the  aforementioned  gamma-rays  (0.05  and  0.1  x rays  per 
alpha  disintegration  respectively  for  J Pu  and  Pu)  (see  Table  i). 


Table  I 

239  240 

Major  Radiation  Energies  and  Yields  in  the  Decay  of  Pu  and  Pu 


Absolute 

240_ 

Absolute 

Radiation 

Enersrv 

Yield 

Radiation 

Yield 

*0 

5.151  MeV 

0.733 

% 

5.168  MeV 

0.76 

*13 

5.1^5 

0.151 

<*45 

5.123 

0.24 

*52 

5.10? 

0.115 

**149 

5.020 

0.0009 

^2 

38.6  keV 

5. 9x1 O'5 

45.23  keV 

4.5X10"4 

K 

51.6 

2.1x10^ 

*2 

103.6 

7.0xl0~5 

^14 

129.27 

6.2xl0-5 

*3 

160.35 

4.2xl0‘6 

1,(0) 

11.618  keV 

0.0011 

Ll(U) 

11.618  keV 

0.0024 

k(u) 

13.615 

0.0182 

ww 

13.615 

0.0378 

i*(u) 

17.220 

0.0216 

i»(u) 

17.220 

0.0484 

Iv(u) 

20.167 

0.0053 

L;(U) 

20.167 

0.0120 

2L(U) 

11.6-20.1 

0.0463 

£L(U) 

11.6-20.1 

0.1006 

6 


Coincidence  Spectroscopy 


As  was  mentioned  in  the  previous  section  on  radioactive  decay 
processes  the  primary  radiation  of  a nuclide  is  often  accompanied  by 
secondary  radiation  of  the  daughter  nuclide.  The  simultaneous  detection 
of  two  or  more  of  these  primary  and  secondary  radiations  by  different 
detectors  can  consequently  yield  valuable  information  about  the  decay 
transitions  and  about  the  excited  energy  states  of  the  daughter  nuclide. 
Coincidence  spectroscopy  is  one  such  method  that  utilizes  the  simul- 
taneity of  detection  to  determine  decay  schemes,  or  decay  rates  of  sam- 
ples. When  used  to  determine  the  decay  rate,  the  method  is  more  appli- 
cable when  applied  to  a radiochemically  pure  sample. 

If  a transition  o(  is  registered  in  one  detector  with  an  efficiency 
, and  in  a second  detector  a transition  J is  registered  with  effi- 
ciency ey  , and  further  assuming  that  both  are  insensitive  to  the  radi- 
ation detected  by  the  other,  then  the  number  of  counts  in  each  detector 
in  a given  time  T for  a sample  with  a constant  decay  rate  R decays  per 
unit  time  is 

N«.-e*ccRT  i (l) 

If  one  further  assumes  isotropic  angular  correlations  between  the 
two  radiations  and  a coincidence  efficiency  of  1.,  the  registered  number 
of  coincidences  recorded  in  the  two  detectors  is 

Nodi'  " ) RT  (2) 

where  <x.  and  j/  are  taken  as  transition  intensities.  This  number 
of  coincidences  includes  coincidences  that  occur  by  change  due  solely  to 
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the  finite  resolving  tine  of  the  coincidence  circuity.  It  can  be  3hovm 
(Ref  3:5*0)  that  the  number  of  random,  or  chance  coincidences  is 

NJa»  E7'(m«.-n«.vX,sJv-n«.v)/t  (3) 

or  when  t\l«.V  is  small  compared  to  N«.  and  hi**  , as  is  often  the  case, 
then 

(4) 

The  true  coincidences  are  essentially  the  observed  coincidences  corrected 
for  chance  coincidence,  dead  time  losses,  and  background. 

The  power  of  the  coincidence  technique  is  that  the  decay  rate  of 
the  sample  can  often  be  calculated  without  direct  knowledge  of  the  in- 
dividual detector  efficiencies.  Combining  equations  (l)  and  (2),  the 
decay  rate,  R,  of  a sample  is  given  by 

R - MotNy/KUj'T  (5) 

The  above  formula  however,  necessitates  that  the  radiations  and  J 
are  both  from  the  same  nuclide.  In  the  case  of  isotopically-mixed  sam- 
ples, this  assumption  is  seldom  true  and  both  transitions  and  J 

are  the  combined  transitions  of  two  or  more  isotopes. 

An  inherent  disadvantage  in  coincidence  spectroscopy  is  that  it  in- 
corporates the  efficiency  of  both  detectors  to  register  a count.  The 
total  efficiency  for  a coincidence  being  the  product  of  individual  de- 
tector efficiencies i consequently,  the  number  of  coincident  counts 
could  be  small  and  the  acquisition  times  large  for  a low-activity  sam- 
ple or  for  a low-geometry  detection  system.  An  advantage  is  that  it 
reduces  background  substantially. 
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239 

In  the  particular  case  of  mixed  samples  containing  both  Pu  and 

240  239 

Pu,  selective  acquisition  of  the  alpha  spectrum  from  Pu  was  at- 
tempted by  gating  on  the  51.6  keV  gamma  ray  emitted  by  this  isotope, 

239 

The  activity  of  ^ Pu  could  then  be  calculated  with  knowledge  of  the 

efficiency  of  the  gamma-ray  detector  and  the  abundance  of  gamma  rays  per 

disintegration.  However,  a problem  arises  in  identification  of  the  gamma 

ray  that  is  to  serve  as  the  pulse  to  trigger  the  gate.  With  abundances 

239 

as  low  as  they  are  for  the  51.6  keV  gamma  ray  of  Pu  it  could  well  be 
undetected  or  hidden  in  the  background. 

Singles  Spectroscopy 

239 

Just  as  the  alpha  spectrum  of  Pu  is  unresolvable  from  that  of 
240 

Pu,  so  also  is  the  uranium  characteristic  x-ray  spectrum  arising  from 

the  decay  of  both  isotopes.  However,  the  fact  that  the  ratio  of  total  L 

240 

x rays  per  alpha  transition  for  Pu  is  approximately  twice  that  ratio 
239 

for  Pu  theoretically  provides  one  with  sufficient  information  from 

the  two  singles  spectra  of  alpha  particles  and  L x rays  to  calculate  the 
239  240 

ratio  of  Pu  to  Pu  in  a given  sample. 

If  one  takes  OCr  as  the  total  alpha  counts  under  the  5.15  MeV  to  5.17 

239 

MeV  alpha  peak  as  being  a sum  of  the  individual  oc  counts  of  Pu  and 
240 

Pu  and  similarly  defines  X/p  as  the  total  L-shell  characteristic  x rays 
from  internal  conversion  of  the  two  isotopes  then 

OCt  * +■  2*0  (6) 

X.T  “ (7) 
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Prom  recent  studies  at  Cak  Ridge  National  Laboratories,  (Ref  4:578), 

and  Argonne  National  Laboratories,  (Ref  5),  the  ratios  of  L x rays  to 

2^9  240 

alpha  disintegration  have  been  well  documented  for  " Pu  and  ?u: 

-**23  =.  0.046  - 0.101  (8) 


By  substituting  these  last  two  ratios  (8)  into  equation  (?)  and 

then  solving  simultaneously  with  equation  (6),  one  obtains  the  activity 
2^q  240 

ratio  of  J Pu  to  Pu  in  terms  of  the  counts  under  the  alpha  peak  and 
under  the  L x-ray  spectrum. 

_ 0,101  <XT  — X-t 

<*.,»  XT—  0.046  t>cT  (9) 


The  atom  ratio  of  these  two  isotopes  can  be  calculated  from  the  activity 
ratio  by  use  of  the  decay  coefficients: 


zw 


Pu 


240 


(10) 
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Ill,  Experimental  Equipment  and  Electronics 


Introduction 

The  initial  thrust  of  this  project  was  directed  to  the  use  of  alpha- 
gamma  coincidence  spectroscopy.  Consequently,  the  selection  of  equip- 
ment and  the  design  of  apparatus  were  influenced  by  the  need  to  max- 
imize efficiencies  and  geometries  for  detection  of  these  particular  ra- 
diations. Because  of  the  low  activity  of  the  samples,  a concentrated 
effort  was  made  from  the  beginning  to  reduce  background  radiation. 

When  later  experimental  efforts  turned  toward  singles  spectroscopy 
of  both  alpha  particles  and  L-series  x rays,  the  equipment  used  for  alpha 
spectroscopy  in  the  original  coincidence  experiment  was  retained  and  a 
detector  and  spectroscopy  circuit  were  incorporated  for  the  x rays. 

Equipment 

The  low  penetration  capabilities  of  alpha  particles  in  air  gener- 
ally requires  that  alpha  spectroscopy  be  accomplished  with  the  sample 
and  detector  in  an  evacuated  chamber.  The  coincident  detection  of  gamma 
rays  from  the  same  sample  requires  that  a second  detector  be  capable  of 
"seeing"  the  sample  simultaneously. 

The  major  design  problem  concerned  the  configuration  of  the  chamber 
and  sample  mount  to  allow  both  alpha  and  gamma  detectors  to  "see"  the 
sample  and  to  maximize  the  solid  angles  for  each  detector. 

Vacuum  Chamber.  Samples  of  heavy  elements  from  McClellan  Central 
Laboratory  (KCL)  consist  of  the  isotopic  material  electrodeposited  on 
3 mil-thick  platinum  disk  of  approximately  2.2  cm  diameter.  The  dia- 
meter of  the  area  of  deposition  is  1.?  cm.  A sample  mount  that  would 
hold  the  disk  by  the  outer  edge  close  to  the  alpha  detector,  and  yet 
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would  not  further  attenuate  the  gamma  radiation  exiting  the  backside  of 

the  sample,  appeared  to  offer  the  best  geometry  for  both  detectors. 

(2 

3ecause  of  the  1/R  dependence  of  solid  angle  on  sanple-to-detector  dis- 
tances, these  distances  were  reduced  to  absolute  minima  allowing  only 
for  slight  variations  in  sample  configuration  and  detector  size.  The 
completed  assembly  consisted  of  an  aluminum  cylinder,  one  inch  high  by 
four  inches  in  diameter.  The  alpha  detector  was  attached  on-axis,  in- 
ternally at  one  end  of  the  chamber.  The  endplate  of  the  chamber  served 
as  sample  mount,  vacuum  seal,  and  exit  window  for  gamma  radiation  and 
had  a beryllium  window  of  7/8  in.  diameter  and  a thickness  of  20  mil. 

A detailed  drawing  of  the  chamber  is  presented  in  Appendix  3, 

With  the  endplate  in  place,  the  active  side  of  the  sample  disk  was’ 
5.5  inm  from  the  face  of  the  alpha  detector.  With  the  gamma  detector 
directly  against  the  chamber,  the  distance  from  sample  to  detector  en- 
trance window  was  4,0  mm.  (The  sensitive  surface  of  the  detector  was 
another  8.0  mm  behind  the  window). 

The  combined  chamber  and  detector  assembly  was  Inserted  within  a 
cylindrical  lead  shield  of  1$  in.  thick  sides,  (see  Fig.  2),  One  end 
of  the  cylindrical  lead  chamber  was  shielded  with  a removable  lead  plate 
of  one  in.  thickness. 

A graded  shielding  consisting  of  a 25-mil-thick  copper  sleeve  was 
inserted  between  lead  shield  and  aluminum  chamber  to  absorb  x rays  in- 
duced in  the  outer  lead  cylinder. 
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Detectors.  The  detector  originally  used  for  alpha  spectroscopy  was 


an  ORTEC,  Model  A-020-100-100  partially-depleted,  silicon  surface- 

2 

barrier  detector  with  an  active  area  of  100  ran  and  a depletion  depth  of 
100/a  . It  was  chosen  for  its  compact  size,  short  charge-collection 
times,  and  excellent  energy  resolution  (total  measured  system  resolution 
was  25  keV  FVHM  at  5*49  KeV),  It  further  provided  for  a geometry  fac- 
tor , fl/qir , of  0.125  (where  H is  the  solid  angle  in  steradians). 

Later,  the  detector  was  replaced  by  a larger  ruggedized  silicon 

surface-barrier  to  increase  the  geometry  factor.  An  ORTEC  Model 

2 

BR-25-^50-100  with  an  active  area  of  450  mm  and  a depletion  depth  of 
100/A.  was  chosen.  Due  partially  to  its  thicker  aluminum  entrance  win- 
dow and  larger  surface  area,  the  system  resolution  was  reduced  to  45  keV 
FVHM  at  5.49  MeV,  The  decrease  in  resolution  did  not  hamper  resolving 
the  alpha  particles  of  interest  and  the  two-fold  increase  in  efficiency 
(geometry  factor  of  0.258)  was  highly  desirable  in  reducing  the  acquisi- 
tion times  for  spectra  of  lower-activity  samples. 

The  detector  selected  for  gamma-ray  spectroscopy  was  a lithium- 
drifted,  germanium,  Ge(Li),  detector  of  right-angle-cryostat-configura- 
tion. The  detector  was  chosen  because  of  its  excellent  energy  and  tim- 
ing resolutions  (system  energy  resolution  was  measured  at  2.8  keV  FVHM 
at  122  keV)  and  further  because  of  its  thin  beryllium  entrance  window, 

providing  for  high  transmittance  of  low-energy  photons.  The  sensitive 

2 

surface  of  the  detector  diode  had  an  area  of  5.0  cm  and  was  mounted 
8,0  mm  behind  the  entrance  window.  With  the  detector  window  physically 
against  the  beryllium  window  of  the  vacuum  chamber,  the  geometry  factor 
for  photon  detection  was  0.139. 
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Electronics 


The  block  diagram  for  the  coincidence  circuit  is  shown  in  Fig.  5 
Appendix  3. 

Amplifiers.  The  preamplifier  for  the  Nuclear  Diode  Ge(Li)  de- 
tector was  integral  to  the  detector  housing  and  provided  connections 
for  detector  bias  voltages  as  well  as  for  signal  inputs  to  the  linear 
amplifier.  The  preamplifier  used  with  the  surface  barrier  detector  was 
an  GRTEC  Model  109A  with  a low-noise  FET  input  stage. 

The  linear  amplifiers  amplify  and  shape  the  output  pulses  of  the 
preamplifiers  to  obtain  optimum  timing  and  energy  resolution.  Tennelec 
Model  TC  203  BLR  linear  amplifiers  were  chosen  for  both  the  alpha  and 
gamma  channels.  This  insured  similar  output  pulse  shapes  and  signal 
processing  times. 

Coincidence  Circuitry.  Bipolar  outputs  of  the  linear  amplifiers 
were  fed  into  timing  single  channel  analyzers  (TSCA)  for  differential 
or  Integral  energy  discrimination  and  also  as  a means  of  inserting 
timing  delays  in  either  channel  prior  to  the  coincidence  unit.  The 
outputs  of  the  TSCA's  were  fed  into  an  ORTEC  Model  4l4A,  fast  coinci- 
dence unit,  which  generates  a logic  pulse  output  if  the  cross  over- 
times of  the  two  input  pulses  are  within  the  resolving  time  of  the  cir- 
cuit (continuously  variable  from  10  nsec  to  110  nsec).  The  output  of  the 
fast  coincidence  unit  was  fed  into  an  ORTEC  Model  41 6A  gate -and -delay 
generator  which  translates  the  logic  pulse  of  the  coincidence  unit  into 
output  pulses  of  selectable  amplitude  and  width.  The  output  pulse  of 
the  gate  and  delay  generator  was  fed  directly  to  the  gate  input  of  the 
multichannel  analyzer  (MCA), 
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Spectroscopy  Circuit 


The  monopolar  output  pulse  of  the  linear  amplifier  for  the  alpha 
channel  passed  through  a delay  amplifier  which  delayed  the  output  suffi- 
ciently to  compensate  for  the  time  required  for  the  coincidence  circuits 
to  process  inputs.  The  output  of  the  delay  amplifier  was  connected  to 
an  ORTEC  Model  408  biased  amplifier  which  expanded  regions  of  interest  in 
the  spectrum  and  its  output  was  routed  directly  into  the  signal  input  of 
the  MCA.  Three  counter  timers  were  used  to  provide  continuous  singles 
counting.  These  were  monitored  to  determine  approximate  sample  activi- 
ties, to  aid  in  trouble  shooting,  and  to  provide  the  singles  rates  for 
the  random  coincidence  corrections  of  equation  (4).  Two  scalers,  RIDL 
Model  49-43,  were  used  to  record  singles  counts  for  the  alpha  and  gamma 
channels.  The  other  counter,  a Canberra  Model  1771,  recorded  total  coin- 
cidences. 

X-Ray  Spectroscopy 

The  spectroscopy  circuit  assembled  for  x rays  is  similar  to  the  cir- 
cuitry for  alpha  spectroscopy,  previously  mentioned,  (see  Fig.  6,  Appen- 
dix B),  The  detector  chosen  was  the  Princeton  Gamma-Tech  Model  LS-33 
lithium-drifted  silicon,  Sl(Li),  detector.  It  was  chosen  because  of  ex- 
cellent resolution  and  high  efficiency  for  low-energy  (<30  keV)  photon 
detection.  The  choice  was  made  between  the  Si(Ll)  and  an  available  in- 
trinsic germanium  detector  primarily  on  the  shorter  distance  from  entrance 
window  to  sensitive  detector  area  for  the  Si(Li)  detector.  This  shorter 
distance  provided  for  a higher  total  efficiency.  Another  factor  in  choos- 
ing the  Si(Li)  over  the  germanium  was  an  observed  lower  background  for 
the  Sl(Li).  Both  detectors  had  active  areas  of  30  mm  and  resolutions  for 
the  two  were  comparables  200  eV  FWHM  3 6.4  keV  for  Ge  and  230  eV  FWHM 
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S 6.4  keV  for  the  Si(Ll)  detector.  The  germanium,  because  of  its  higher 
atomic  number,  provided  for  higher  detection  efficiencies  for  photons 
greater  than  30  keV,  where  the  efficiency  of  Si(Li)  drops  off  rapidly, 
but  since  the  L x rays  of  uranium  range  in  energy  from  11-21  keV,  this 
was  of  no  concern. 

The  preamplifier  of  the  x-ray  detector  was  built  into  the  cryostat 
assembly  and  its  output  was  fed  to  a Princeton  Gamma-Tech  Model  340 
linear  amplifier.  The  output  of  the  linear  amplifier  was  fed  into  a 
Tennelec  Model  TC601  strobed  pulse  stretcher,  which  provided  for  better 
energy  resolution  in  shaping  the  pulse  to  that  compatible  with  the  MCA. 

A biased  amplifier  was  inserted  between  pulse  stretcher  and  MCA  to  expand 
the  11  keV  to  21  keV  spectrum  of  interest  to  a full  scale  range  on  the 
MCA. 


Figure  3.  Sample  Mount  for  X-Ray  Detector 
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The  x-ray  detection  system  required  a sample  mount  that  would  hold 
the  source  against  the  detector  housing  in  a reproducible  geometry.  The 
final  design  consisted  simply  of  a nylon  sleeve  that  fit  snugly  over  the 
detector  housing.  The  sample  was  held  onto  the  end  of  the  sleeve  by  a 
l/l6  in.  thick  piece  of  lead  wrapped  around  the  outside  of  the  sleeve. 

The  front  of  the  detector  was  shielded  by  an  additional  piece  of  l/8  in. 
thick  lead.  The  entire  assembly  was  of  light  weight,  provided  sufficient 
shielding  for  low-energy  photons  and  yet  permitted  an  easy  change  of  sam- 
ples. (see  Fig.  3) 
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IV,  Experimental  Procedures  : Calibration 
and  Data  Collection 

Several  calibration  measurements  were  necessary  before  raw  spectral 
data  could  be  identified  or  quantified,  Kost  calibrations  were  a one- 
time accomplishment,  still  others  had  to  be  checked  periodically  for 
changes  in  system  response.  Whenever  any  change  in  a component  was 
deemed  necessary,  a careful  analysis  of  the  impact  on  the  system  was 
essential  to  assess  the  necessity  of  a recalibration. 

Energy  Calibration 

Calibrations  are  necessary  to  relate  energy  to  channel  number  for 
the  pulse  height  spectra  collected  by  the  MCA.  This  is  generally  a one- 
time measurement  for  each  detector-amplifier  combination}  however,  it 
must  be  repeated  with  any  change  in  detector,  amplifier  component,  de- 
tector bias,  or  amplifier  gain.  The  energy  calibration  techniques  of 
the  system  for  alpha,  gamma,  and  x-ray  radiation  are  all  quite  similar, 
in  that,  all  use  sources  emitting  radiation  of  referenced  energies.  The 
choice  of  particular  sources  used  were  dependent  upon  the  full-scale 
energy  range  of  interest. 

Alpha  Spectroscopy.  Several  alpha  reference  sources  were  available 
o/li  oqq  T.'VK 

for  energy  calibration:  iAm,  J U,  however,  the  most  useful 

source  for  much  of  this  work  was  an  isotopically  mixed  sample  from  KCL 
24?  243  244 

containing  j?u,  Am,  and  Cm  and  providing  alpha  energy  peaks  of 
about  4.9,  5.3,  »nd  MeV  respectively.  This  sample  contained  several 
thousand  DPM  of  each  nuclide  electrodeposited  on  the  standard  Pt  disk  and 
allowed  rapid,  high-resolution  determination  of  peak  channels.  These 
reference  peaks  spanned  the  energy  range  of  interest  very  well. 
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It  was  possible,  after  initial  energy  calibration,  to  maintain  a 

continuous  cross  check  of  energy  calibration  by  reference  to  the  peaks 
239  240  241 

of  Pu,  Pu  and  Am  normally  present  in  each  sample  being  analyzed. 

Gamma  Spectroscopy . For  gamma-ray  energy  calibration  the  range  of 

energies  being  analyzed  was  generally  confined  to  that  less  than  150 

57  109 

keV,  The  most  convenient  sources  were  found  to  be-  those  of  J Co,  Cd, 
241 

and  Am,  These  three  isotopes  provide  photon  reference  peaks  of  14.4, 
22,  59.5,  87.7,  122  and  136  keV. 

X-Ray  Spectroscopy.  For  initial  energy  calibration,  the  sources  of 
57  10° 

Co  and  'Cd  providing  photon  energies  of  6.4,  7.05,  14.4,  and  22  keV 

from  gamma  transitions  and  characteristic  x rays  were  most  useful. 

After  expansion  of  the  spectrum  to  the  energy  range  of  Interest  1 (11-21 

keV)  by  the  biased  amplifier,  the  L«_  , Lyj  , and  Ly  x rays  from  internal 

24l 

conversion  in  Np  following  the  decay  of  Am  were  used  for  follow-on 
monitoring  of  energy  calibration. 

To  these  several  reference  points  of  energy-vs. -channel-number  for 
each  of  the  above  three  spectra,  a least-squares,  straight  line  was  fit 
using  the  linear  regression  routine  of  a Texas  Instrument  SR-51  hand 
calculator.  The  slope  of  the  line  and  intercept  were  used  to  determine 
the  energy  of  peaks  so  they  could  later  be  identified. 

Efficiency  Calibration 

The  method  for  efficiency  calibration  of  alpha,  gamma  and  x-ray  de- 
tectors is  essentially  the  same,  in  that  all  were  calibrated  by  reference 
to  samples  of  a known  activity.  Each  detector  calibration,  however, 
offered  its  own  particular  problem  and  difficulty. 
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Alpha  Particle  System.  The  large  ruggedized  surface-barrier  detec- 


tor was  calibrated  for  efficiency  by  reference  to  the  supplied  activity 
figures  of  three  separate  samples  supplied  by  MCL.  These  were  mounted 
in  the  standard  geometry  configuration  in  the  chamber.  Table  1 lists  the 
sample  by  number  and  the  activity  for  each  isotope  present.  The  uncer- 
tainties in  the  calculated  efficiencies  are  the  standard  deviations  and 
include  the  uncertainty  in  sample  activity  as  well  as  the  statistical  un- 
certainty in  the  measured  count  rate. 

Table  II 

Calculated  Total  Efficiencies  for  Standard  Alpha  Sources 


Sample  # 

Nuclide 

Alpha  Enerjrv  (MeV) 

Corrected 
Activity  (DPM) 

Total 

Efficiency 

Cm-1 

242Pu 

4.90 

3214+15 

0.255+. 002 

243A* 

5.27 

3182+15 

0.257+. 003 

5.81 

1894+11 

0.249+.003 

Pu-02 

242Pu 

4.90 

632+6 

0.270+.003 

239pu 

5.15 

689+7 

0.263+.003 

5.49 

54+2 

0.283+.010 

Pu-03 

239pu 

5.15 

460+6 

0. 254+.  002 

2*Pu 

5.77 

405+5 

0. 277+. 003 

The  above  efficiencies  appeared  to  have  no  energy  dependence  and  were 
consequently  averaged  to  give  a value  oft  0. 264+. 004.  The  intrinsic  effi- 
ciency for  alpha  particles  incident  on  the  face  of  the  detector  was  ini- 
tially assumed  to  be  very  close  to  100^. 
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A program  SCLIDA  was  written  that  would  calculate  the  solid  angle 
for  a circular  disk  source  on  axis  and  at  a given  distance  from  a cir- 
cular disk  detector,  (see  Appendix  C).  This  geometry  factor  was  then 
compared  to  the  measured  efficiency.  Using  the  area  of  sensitive  region 
of  the  detector,  the  area  of  deposition,  and  the  source-to-detector  dis- 
tance! the  geometry  factor  was  calculated  to  be  0.258+. 008.  This  figure 
was  only  slightly  lower  than  the  measured  total  efficiency  of  0.26*4-+. 004. 
The  difference  was  within  expected  statistical  variations)  however,  the 
factor  of  1.02  difference  in  the  two  was  examined  for  a backscattered 
contribution  from  the  platinum  backing.  A.  H,  Snell  gives  values  as 
high  as  1.04  for  backscattered  alphas  from  a Pt  backing  when  detected  in 
a 2lf  geometry  (Ref  13:336).  Siegbahn,  however,  states  the  contribution 
from  backscatter  in  a point  source-detector  geometry  of  steradians  is 
negligible,  since  most  backscattered  alpha  particles  depart  the  surface 
at  low  angles  (Ref  3:428).  Since  the  sample-detector  geometry  for  this 
experiment  cannot  be  represented  by  a point  source,  some  backscattered 
alpha  particles  may  be  incident  on  the  detector. 

Gamma-Ray  System.  The  problem  of  efficiency  calibration  for  photon 
detection  is  not  as  easily  handled  as  that  for  alpha  particle  detection. 
The  problem  arises  from  the  highly  dependent  nature  of  photoelectric  ao- 
eorption  on  photon  energy.  As  the  energy  of  the  photon  increases,  the 
photoelectric  coefficient  decreases.  On  the  other  hand,  as  the  photon 
energy  decreases,  the  more  it  is  attenuated  by  the  detector  window,  air 
between  sample  and  detector,  and  by  self  absorption  of  the  sample  mate- 
rial itself.  Consequently,  photon  total  efficiency  curves  are  charac- 
terized by  curves  that  drop  off  at  both  ends  of  the  energy  spectrum. 
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For  the  particular  case  of  the  gamma  rays  In  this  experiment,  the 
area  of  interest  was  a narrow  range  of  energies  on  either  side  of  the 
51.6  keV  gamma  ray  of  Pu.  This  range  was  primarily  limited  to  40  to 
150  keV  and  gamma-ray  standards  supplied  by  N3S  were  sufficient  to  de- 
fine the  calibration  curve  over  this  range.  The  photopeak  efficiency 
57  109  241  137 

curve  derived  from  Co,  Cd,  Am,  and  Cs  after  corrections  were 
made  for  source-geometry  differences  and  window  absorption  is  plotted  in 
Fig.  7. 

X-Ray  System.  Accurate  efficiency  calibration  for  the  x-ray  detec- 
tor was  most  difficult.  Much  of  the  uncertainty  in  the  final  results  of 
this  experiment  is  due  to  the  uncertainty  of  this  efficiency.  The  prob- 
lem primarily  arises  as  a result  of  the  scarcity  of  standardized  samples 

which  emit  gamma  rays  or  x rays  in  the  lower  energy  ranges  (<40  keV). 

57  1 37  1 09 

The  gamma -ray  standards,  J Co,  Cs,  and  Cd,  are  primarily  used  to 
calibrate  gamma-ray  detectors  over  a higher  energy  of  50-660  keVj  how- 
ever, these  sources  also  emit  x rays,  the  intensities  of  which  have  to 
be  derived  from  internal  conversion  coefficients  and  fluorescent  yields. 
Larger  uncertainties  have  to  be  attached  to  these  emission  rates  than 
that  of  the  primary  gamma  rays.  Moreover,  since  the  sources  are  prim- 
arily higher  energy  gamma  standards,  the  sandwiching  material  of  the 
sample  can  greatly  effect  the  emission  rates  of  the  lower  energy  photons. 
(Ref  6t20?-208),  The  small  entrance  window  of  the  Sl(Li)  detector  made 
sample  alignment  very  critical.  The  areas  of  deposition  for  these  gamma- 
ray  standards  were  nominally  2-3mm  in  diameter  but  small  shifts  off  axis 
could  greatly  effect  the  geometry  and  resulting  efficiency  data. 

The  most  accurate  approach  appeared  to  be  that  of  measuring  charac- 
teristic x rays  of  Np  from  internal  conversion  following  the  decay  of 
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(Ref  6:215).  Two  americium-241  standards  issued  by  N3S  and  NEN 
(New  England  Nuclear  Co.)  were  employed  for  this  measurement.  These 
sources  were  better  than  the  aforementioned  gamma-ray  sources  because 
as  alpha  sources  they  were  not  sandwiched  in  polyethylene.  The  stan- 
dardized activities  of  the  samples  were  decay  corrected  and  the  abso- 
lute intensities  of  x rays  per  disintegration  used  were  those  of 
Campbell  and  McNelles  (Ref  6:215). 

After  correcting  for  geometry  differences  and  window  absorption, 
the  resulting  efficiencies  were  plotted  vs,  energy  with  those  determined 
by  the  gamma-ray  standards  (see  Fig.  8).  The  discrepancy  in  the  two 
efficiency  curves  derived  from  the  two  types  of  standards  were  never 
resolved  despite  repeated  measurements. 

The  efficiency  data  derived  from  the  Np  x rays  were  used  to  quan- 

239  240 

tify  x ray  abundances  of  Pu  and  Pu  in  this  experiment.  It  was 
felt  that  the  very  small  energy  difference  in  characteristic  x rays  of 
Np  and  U v>ould  introduce  less  error  for  efficiency  than  that  of  the 
Interpolated  curve  derived  from  the  wider  spread  energies  of  the  gamma 
ray  standards. 

Relative  Delay  Calibration 

The  relative  delay  of  the  gamma-ray  and  alpha  particle  detector 
signals  being  fed  to  the  coincidence  unit  is  of  prime  importance  in  de- 
termining the  total  coincidence  count  rate.  These  two  signals  must 
arrive  at  the  coincidence  unit  within  a very  short  time  interval  to  in- 
sure operation  at  the  maximum  coincidence  rate.  Failure  to  operate  at  or 
near  this  maximum  rate  significantly  decreases  the  total  number  of  coin- 
cidences and  increases  the  relative  contribution  of  random  coincidences 
to  the  total  coincidence  count. 
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This  optimum  relative  delay  was  determined  by  varying  the  delay  of 
the  signal  in  one  branch  of  the  coincident  circuit  by  the  delay  poten- 
tiometer of  the  TSCA.  A source  that  emitted  a relatively  high  number 
of  alpha  particles  in  coincidence  with  gamma  rays  was  placed  in  the 
vacuum  chamber  for  simultaneous  alpha  particle  and  gamma-ray  detection. 

A curve  of  coincident  counts  vs.  relative  delay  was  then  plotted  for  a - 
110  nsec  resolving  time  setting  on  the  fast  coincidence  unit,  (see 
Fig.  9)*  The  measure  of  FWHM  of  the  delay  curve  provides  the  resolving 
time,  2T,  of  the  coincidence  circuit.  It  is  advantageous  to  make  the 
resolving  time  as  small  as  possible  to  reduce  the  number  of  chance  coin- 
cidences, It  should  not,  however,  be  reduced  to  a value  that  would  lower 
the  coincidence  efficiency  less  than  100 When  the  delay  curve  has  a 
flat  top,  one  can  usually  feel  confident  that  the  coincidence  efficiency 
will  be  nearly  equal  to  one,  (Ref  3 *55*0 » 

Data  Collection 

For  coincidence  spectroscopy  when  gating  off  a particular  energy 
gamma  ray,  the  upper  and  lower  level  discriminators  of  the  TSCA  must 
first  be  set  to  estaolish  the  energy  window.  If  the  window  is  set  too 
narrow,  true  coincidences  will  be  lost.  If  the  window  is  too  wide  or 
displaced  from  the  gamma -ray  energy  of  interest,  an  excessive  number  of 
unwanted  coincidences  will  be  counted.  This  window  is  set  by  placing  the 
TSCA  in  differential  mode  and  routing  its  output  directly  to  the  gate  in- 
put of  the  MCA.  This  permits  gating  on  the  same  spectrum  being  acquired 
by  the  MCA.  By  varying  the  lower  level  discriminator  and  AS  during  ac- 
quisition of  a gamma-ray  spectrum,  the  spectrum  can  be  limited  to  that 
of  the  desired  energy  window. 
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After  setting  the  energy  window,  the  acquisition  of  the  gated  spec- 
trum is  simple  and  only  entails  setting  of  the  desired  acquisition  time. 
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V. 


Data  Reduction  and  Srror  Analysis 


Coincidence  Sue ctros cony 
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The  very  low  Intensity  of  gamma  rays  emitted  in  the  decay  of  Pu 
240 

and  Pu,  plus  the  low  activity  levels  of  the  samples  themselves,  con- 
tributed to  the  unfeasibility  of  alpha-gamma  coincidence  spectroscopy 
at  these  activity  levels.  The  number  of  coincident  counts  recorded  in 
any  one  channel  over  a 48-hour  period  was  never  more  than  three  >•  conse- 
quently, there  were  not  any  significant  data  to  be  reduced.  Data  re- 
duction for  these  spectra  consisted  of  calculating  the  contribution  of 
random  coincidences  to  the  total  coincidence  spectrum.  This  was  done 
by  using  equation  (4),  Chapter  II,  to  calculate  the  predicted  number  of 
coincident  counts  in  the  acquisition  period.  Had  this  number  ever  been 
significant,  a random  coincidence  spectrum  would  have  been  obtained. 

This  would  have  been  accomplished  by  setting  the  relative  delay  in  one 
branch  of  the  circuitry  approximately  .5  ><sec  off  of  the  delay  curve 
plateau.  The  acquired  spectrum  would  then  have  been  due  to  chance  co- 
incidences only. 

Singles  Sue ctros copy 

Quantification  of  a given  radiation  emitted  by  a sample  consisted 
of i totaling  the  number  of  counts  under  the  energy  peakj  subtracting  the 
background  contribution j then  dividing  this  net  count  by  the  total  effi- 
ciency of  the  detector. 

Alpha  Spectra.  Determining  the  number  of  alpha  particles  emitted 
by  a sample  in  a given  period  of  time  was  a simple  process  of  subtract- 
ing a pre-measured  background  spectrum  from  the  gross  count  spectrum  and 
dividing  by  the  source  detector  efficiency.  The  resolution  of  the 
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detector,  and  the  quality  and  uniformity  of  the  samples  were  sufficient 

2^1  238 

that  the  energy  peaks  of  am  and  Pu  were  resolvable  from  the  com- 
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bined  J^?u/  **  Pu  energy  peak.  It  was,  therefore,  unnecessary  to  cor- 
rect for  tailing  contribution  to  adjacent  peaks  as  in  the  case  of  low- 
resolution  systems  (Ref  9 * 1352) . The  background  had  to  be  periodically 
measured  throughout  the  experiment  because  of  increases  in  chamber  and 
detector  contamination  from  the  Pu  samples. 

X-Ray  Spectra.  Determination  of  the  net  number  of  counts  under  an 
energy  peak  of  the  x-ray  spectrum  was  more  involved  than  that  for  alpha 
spectroscopy.  The  background  for  photon  spectroscopy,  is  in  general, 
very  source  dependent.  It  includes,  in  addition  to  natural  background, 
that  due  to  scattered  radiation  of  the  source  Itself, 

There  was  neither  a simple  not  accurate  method  of  predetermining 
the  background  spectrum  prior  to  the  actual  acquisition  of  the  sample 
spectrum.  The  background  was,  instead,  determined  in  a paired  observa- 
tion with  the  sample  spectrum.  For  the  case  of  the  lower  activity  sam- 
ples (''■'1000  DPM),  the  acquired  peak  counts,  as  well  as  background,  were 
low  for  even  extensive  counting  periods  of  48  hours  (see  Fig,  11).  The 
background  counts  in  these  cases  were  averaged  over  a twenty- channel  span 
at  five  or  six  places  throughout  the  spectrum  between  peaks.  It  was  de- 
termined that  this  background  could  be  approximated  very  well  by  a 
straight  line,  A least-squares  straight  line  was  fit  to  these  five  or 
8 lx  background  points  and  the  line  was  subtracted  from  the  gross  spectrum 
producing  the  net  number  of  counts  under  each  peak.  Poisson  statistics 
were  used  in  calculating  the  uncertainties  in  background  and  peak  counts, 
and  the  error  was  propagated  in  subtraction.  (Ref  12 156), 
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The  above  procedure  was  satisfactory  for  determining  the  number  of 

counts  under  a U L x-ray  peak  for  an  elementally  pure  sample.  Often, 

however,  the  spectrum  consisted  of  unresolved  complex  peaks  from  the 
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contribution  of  other  actinide  characteristic  x rays.  Am  is  one  such 

24 1 

isotope  often  found  in  Pu  samples.  The  alpha  decay  of  Am  is  followed 

by  the  emission  of  Np  characteristic  x rays.  The  L x ray  intensity  for 

241  240 

Am  is  nearly  four  times  that  of  Pu  and  more  than  eight  times  that 

of  2^pu  (Rex'  4:530),  If  the  difference  in  decay  rates  of  these  isotopes 

is  considered,  the  contribution  of  Np  x rays  to  the  compound  Np/u  x-ray 
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spectrum  could  be  significant  with  only  a small  atom  abundance  of  Am 

in  the  sample.  Alpha  particle  spectroscopy  can  be  used  to  determine 

2i+*Am  content  unless  23®Pu  is  present  too.  The  alpha  groups  of  these  two 

isotopes  are  also  unresolvable , In  the  past,  gamma-ray  spectroscopy  has 
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been  vised  effectively  in  determining  Am  content  in  a Pu  sample 

(Ref  7 *198)  (Ref  8:1759) l however  it  has  necessitated  the  calibration  of 

a detector  for  gamma  rays.  The  detection  efficiency  for  59  keV  photons 

by  a Si(Li)  detector  is  very  low,  so  this  would  effectively  require  a Ge 

or  Ge(Li)  detector.  A method  for  spectrally  subtracting  the  NP  L x-ray 

contribution  from  the  compound  spectrum  would  be  valuable  in  accurately 

quantifying  U L x rays  when  only  a Si(Li)  detector  is  in  use. 

In  spite  of  the  compound  nature  of  the  NP/U  L x-ray  spectrum,  the 

lA»  peaks  of  NP  and  U were  sufficiently  resolved  that  a spectrum  of 

gaussian  peaks  could  be  fit  to  the  individual  Np  L x-ray  lines  based  upon 

the  intensity  of  the  Np  1^,  peak.  This  theoretical  spectrum  for  Np  L 
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x rays  following  Am  decay  was  generated  by  program  SPECTRA  ( see 
Appendix  D)  using  energies  and  relative  intensities  of  Np  L x-ray  sub- 
group lines  provided  by  Hyde  (Ref  11:253*254),  This  same  program 
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performs  a third -degree , least-squares  polynomial  fit  to  the  background, 
adds  the  background  and  Np  L x-ray  spectrum  then  substracts  that  sum 
spectrum  from  the  compound  Np/U  L x-ray  spectrum.  The  remaining  spectrum 
is  that  due  only  to  U L x rays.  The  counts  under  these  peaks  were  then 
totaled  and  quantified. 

The  above  procedure  of  spectra  "stripping"  was  necessitated  by  the 
lack  of  literature  on  individual  U L x-ray  line  intensities  following  Pu 
alpha  decay.  The  subgroup  intensities  of  L*  t L*.  , L^g  and  Lj)  have 
been  documented  (P.ef,  5)  (Ref  4:580)  and  the  energies  are  well  known 
but  subgroup  x-ray  line  intensities  are  needed  to  predict  a theoretical 
U x-ray  spectrum. 

CURVFIT,  an  interactive  graphics  program  that  fits  either  linear  or 
nonlinear  functions  to  a set  of  data  points  was  used  in  an  attempt  to  de- 
termine its  utility  for  quantifying  U L x rays  in  a complex  spectrum. 

Its  major  limitations  were  that  it  could  handle  only  200  channels  of 
data  points  at  one  time  and  was  limited  to  a function  of  ten  or  fewer 
independent  variables.  With  energies  and  standard  deviations  specified, 
this  would  allow  ten  Gaussian  peaks  to  be  fit  by  amplitude  variations  or 
only  five  peaks  if  both  amplitude  and  energy  were  varied  in  the  fit. 

With  the  amplitudes  supplied  by  the  program,  the  areas  under  the  peaks 
could  then  be  calculated  by  a separate  program.  Because  of  the  1024- 
channels  used  for  the  x-ray  spectra  and  the  complexity  of  peaks  that 
must  be  fit,  the  program  was  found  to  be  only  marginally  useful. 

Error  Analysis 

The  uncertainties  associated  with  quantification  of  radiation  depend 
upon  the  statistical  uncertainties  of  the  counting  process,  as  well  as 
those  of  source-detector  efficiency.  The  first  of  these  uncertainties 
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has  been  briefly  covered  in  the  preceding  section*  the  second  was  more 
difficult  to  determine. 

The  geometry  factor  as  calculated  by  SOLIDA  depends  upon  three  meas- 
urements: source  radius,  detector  radius,  and  source-to-detector  dis- 
tance, The  radius  of  the  detector  was  calculated  from  the  sensitive 
area  of  detector  supplied  by  factory  specifications.  The  radius  of  sam- 
ple deposition  was  supplied,  in  the  case  of  KCL  samples,  and  measured 
for  others.  The  source-to-detector  distances  were  measured  with  a ruler. 
The  95^  confidence  interval  of  certainty  was  assumed  to  be  within  +0,5mm. 
The  uncertainties  associated  with  detector  and  sample  radii,  when  sup- 
plied, were  assumed  to  be  small  in  comparison  with  the  uncertainty  of 
source-to-detector  distance. 

Geometry  factors  were  calculated  by  SOLIDA  for  both  the  measured 
distances  and  for  +C.5mm  of  the  measured  distances.  The  resulting  range 
of  efficiencies  were  taken  as  +2c'  of  the  actual  efficiency.  This  un- 
certainty was  then  propagated  throughout  future  calculations  involving 
the  geometry  factor. 

Minimum  Levels  for  Detection  and  Quantification 

The  usefulness  of  a procedure  for  the  analysis  of  radioactivity  is 
dependent  upon  the  level  of  activity  to  which  it  is  sensitive.  It  was 
felt  that  the  results  of  this  experiment  would  be  more  pertinent  if  min- 
imum activity  levels  for  detection  and  quantification  were  supplied  for 
each  method. 

Numerous  approaches,  some  of  which  are  Inconsistent,  that  define 
these  limits  are  described  in  the  literature.  Currie  of  the  National 
Bureau  of  Standards  (Ref  14*586)  provides  a thorough  discussion  and  work- 
ing formulae  for  calculating  these  limits.  Briefly,  he  defines  three 
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specific  levels  of  interest  (l)  Lc,  the  "decision  limit"  at  which  one 
may  decide  whether  or  not  the  results  of  the  analysis  indicates  detec- 


tion, (2)  Lj,  the  "detection  limit"  at  which  a given  analytical  pro- 
cedure can  be  relied  upon  to  lead  to  detection,  and  (3)  L^,  the  "deter- 
mination limit"  at  which  a given  procedure  will  be  sufficiently  precise 


to  yield  a satisfactory  quantitative  estimate.  The  degree  to  which  the 
detection  limit  can  be  relied  upon  is  established  by  a user-supplied 
acceptable  confidence  interval.  The  precision  of  the  quantitative  re- 
sults determined  at  is  similarly  defined  by  the  standard  deviation  the 
user  is  willing  to  accept  in  the  results. 


Table  III 

Minimum  Levels  for  Detection  and  Quantification 

Lc  (counts)  (counts)  Lq  (counts) 

Paired  Observation  2.33 J/tT  2.71+4.6 5/^  5o(l+fl+ 

(<*:->«.)  1 L ,WJ  j 

"Well-known"  Blank  1.6 2.71+3.29v^  50<l+[l+ 

(&.•  0 ) v J 


These  limits  are  naturally  dependent  upon  the  background  and  its  un- 
certainty. Currie  defines  the  "blank"  as  the  signal  resulting  from  a 
sample  which  is  identical,  in  principle,  to  the  sample  of  interest,  ex- 
cept that  the  substance  sought  is  absent  (or  small  compared  to  the  stan- 
dard deviation  of  the  blank).  The  blank  as  defined  here  is  the  back- 
ground and  represented  by  the  true  mean  of  the  blank)  and  its  stan- 
dard deviation  <5S  . If  the  blank  is  determined  in  a paired  observation 
with  the  signal,  <5sis  not  small  in  comparison  to  the  deviation  of  the 
gross  (signal  + background)  signal rr<Cll. 
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On  the  other  hand,  if  the  background  has  been  determined  in  a separate 
measurement  to  a high  degree  of  precision,  then  the  blank  may  be  con- 
sidered "well-known"  and  is  small  in  comparison  to  <5^.  For  radio 

activity,  observations  are  assumed  to  be  governed  by  the  Poisson  dis- 
tribution, Table  III  provides  the  working  expressions  for  Lc,  L^,  and 

L for  radioactive  analysis.  The  confidence  is  95f»  in  making  the  de- 

<1 

cision,  "detected  or  not  detected"  and  the  relative  standard  deviation 
in  the  quantitative  results  at  the  level  is  10/5. 
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VI.  Results.  Conclusions  and  Recommendations 


The  goal  of  this  work  was  to  explore  methods  and  examine  feasi- 
bilities of  using  high-resolutlon  alpha  and  photon  spectroscopy  to  re- 
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solve  the  two  isotopes,  Pu  and  Pu.  The  thrust  of  the  effort  was 
initially  toward  that  of  coincidence  spectroscopy,  however,  the  empha- 
sis was  later  shifted  to  single  alpha  particle  and  x-ray  spectroscopy. 
This  goal  has  been  partially  obtained.  Alpha-gamma  coincidence  spectro- 
scopy has  been  determined  unfeasible  at  the  activity  levels  proposed  in 
the  scope  of  this  experiment.  Minimum  levels  for  detection  and  quanti- 
fication have  been  theoretically  established  for  the  coincidence  tech- 
nique when  using  the  equipment  assembled  for  this  experiment.  The  re- 
sults of  alpha  and  L-shell  x-ray  singles  spectroscopy  are  promising. 

Due  to  a poor  sample-detector  geometry,  the  efficiency  calculations  for 
x-ray  detection  have  large  associated  uncertainties.  These  uncertainties 
plus  statistical  uncertainties  from  the  small  number  of  counts  have  con- 
sistently produced  L x-ray  to  alpha  disintegration  ratios  with  errors 
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too  high  for  quantitative  analysis  of  Pu/  Pu  ratios.  Due  to  time 
limitations  of  the  experiment,  these  uncertainties  have  not  been  reduced 
and  possible  sources  of  other  systematic  errors  have  not  been  dis- 
covered. A valuable  and  large  contribution  to  this  analysis  has  been  the 
approach  and  methodolology  established  in  this  experiment. 

Coincidence  Spectroscopy  Discussion 

To  provide  the  highest  geometry  for  both  alpha  and  gamma-ray  detec- 
tion, the  gamma  rays  were  detected  through  the  backside  of  the  sample. 
Even  the  3 mil- thick  Pt  disks,  on  which  MCL  routinely  electrodeposit 
their  samples,  attenuate  40-50  keV  photons  greatly. 
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The  transmittance  is  only  0.485  and  0.339  for  52  and  45  keV  photons  re- 
spectively. An  improvement  in  transmittance  however,  was  provided  by  a 
switch  in  backing  materials  to  that  of  5 mil-thick  stainless-steel.  The 
two  transmittances  were  increased  to  0,84  and  0,76  for  52  and  45  keV 
photons  respectively. 

Still  higher  transmittances  were  provided  by  samples  electro- 
deposited  on  metallic  foils  and  metallically-coated  mylar.  Sample  back- 
ings of  gold  and  nickel  coated  mylar,  and  foils  of  stainless-steel  and 
copper  were  used  to  assess  the  merits  of  each.  The  Au-coated  mylar  had 
problems  of  flaking  and  uneven  deposition,  however,  the  Ni-coated  mylar 
accepted  the  deposition  quite  evenly  with  very  little  noticeable  ten- 
dency toward  flaking.  The  stainless-steel  and  copper  foils  of  0,l6  mil 
2 2 

(3.19  mg/cm  ) and  0.15  mil  (3.41  mg/cm  ) thickness  respectively  provided 
high-quality,  evenly-deposited  samples.  Based  upon  visible  quality  and 
transmittance,  the  slight  edge  would  go  to  the  use  of  stainless-steel 
foil.  The  transmittance  for  50  keV  photons  was  0.994  and  that  for  45  keV 
photons  was  0.989.  The  copper  foil  provided  for  transmittances  of  0.993 
and  0.988  respectively  for  52  and  45  keV  photons. 
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Gamma-ray  spectra  were  taken  of  the  front  sides  of  ^ Pu  and  Pu 

samples  in  an  attempt  to  identify  the  gamma  ray  that  was  to  serve  as  the 

gating  pulse  for  the  alpha  spectra.  Repeated  measurements  failed  to 

identify  the  gamma  ray  of  Interest.  The  USAF  Radiological  Health  Lab 
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ran  one  ^ Pu  sample  of  approximately  500  DPM  overnight  on  their  planar 
Ge(Li)  detector  in  a very  low  background  chamber  and  were  unable  to  iden- 
tify the  52  keV  gamma  ray.  Various  attempts  at  gating  off  of  the  area 
of  the  spectrum  corresponding  to  52  keV  yielded  no  significant  coinci- 
dent counts. 
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The  reverse  coincidence  experiment  of  gating  off  the  5tl5  KeV  alpha 
group,  when  acquiring  the  gamma-ray  spectrum  of  the  sample,  was  per- 
formed in  an  attempt  to  reduce  the  background  sufficiently  to  see  the 
gamma-ray  peak.  These  coincidence  experiments  permitted  identification 
of  the  peak,  but  with  only  3 counts  under  the  peak,  it  could  not  be 
quantitatively  analyzed. 

Higher  activity  levels  (1500  DPM)  naturally  produced  more  coin- 
cidences, but  in  a 48  hour  acquisition  period  the  number  of  coincident 
counts  under  a peak  were  still  far  too  few  to  quantify  the  isotopic 
ratio. 

Coincidence  Spectroscopy  Results 

Based  upon  the  total  efficiencies  for  alpha  particle  and  gamma-ray 

239 

detectors  used  in  this  experiment,  the  gamma-ray  intensities  of  Pu 
240 

and  Pu,  and  upon  the  background}  the  minimum  activity  levels  for  de- 
tection and  quantification  by  this  coincidence  technique  can  be  theo- 
retically calculated.  The  average  background  count  ) for  coinci- 
dence spectroscopy  over  even  extended  counting  periods  of  48  hours  was 
zero.  Assuming  the  sample  is  deposited  on  one  of  the  thin  foils  pro- 
viding approximately  99%  transmittance  for  52  and  45  keV  photons,  the 

activity  necessary  (in  DPM)  for  detection  and  quantification  of  a pure 
239  240 

Pu  or  Pu  source  are  tabulated  in  Table  IV  for  various  acquisition 
times. 

Singles  Spectroscopy  Discussion 

Prom  equation  (9),  Chapter  II,  one  observes  that  the  accuracy  of 
this  method  is  very  sensitive  to  the  calculated  x-ray  to  alpha  emission 
ratio  of  the  sample.  Dividing  both  numerator  and  denominator  of  the 
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right  hand  side  of  this  equation  by  <XT,  the  equation  directly  incor- 
porates these  ratios: 


PC  hi 


°-101  --qc; 
~k;  “°‘ou6 


(11) 


This  ratio  of  total  L x-rays  to  alpha  particles  emitted  will  vary  theo- 

, 239  240 

retically  from  0.046  for  pure  Pu  to  0,101  for  pure  Pu.  A calculated 

ratio  outside  this  range  gives  meaningless  results.  This  experiment  has 

been  essentially  concerned  with  measuring  these  ratios  for  the  three 

major  L-shell  x-ray  subgroups:  1^  , , L,>  . The  samples  measured, 

varied  from  400  DPM  to  1600  DPM.  These  samples,  in  some  cases,  pos- 

24i  238 

sessed  carrier  isotopes  and  impurities  of  either  iAm  or  J Pu,  or 
both.  Corrections  to  the  L x-ray  yields  for  their  presence  introduced 
still  other  errors  and  uncertainties  in  the  final  results. 


Table  IV 

Minimum  Sample  Activities  for  Detection  and 
Quantification  in  Coincidence  Spectroscopy 

(2,71  Counts)  (100  Counts) 


Acquisition 
Time  (Hrs) 

239Pu  (DPM) 

240Pu  (DPM) 

239pu  (DPM) 

240Pu  (DPM) 

1 

10,893 

4,728 

401,977 

174,480 

4 

2,723 

1,182 

100,494 

43,620 

8 

1,362 

592 

50,247 

21,810 

16 

681 

295 

25,123 

10,905 

24 

453 

197 

16,749 

7,270 

The  accuracy  of  the  measured  L x-ray  to  alpha  ratios  were  depen- 
dent upon  the  total  efficiencies  and  associated  uncertainties  for  the 
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alpha  particle  and  x-ray  sample-detector  combinations.  Because  of  the 
fact  that  the  area  of  deposition  was  larger  than  the  SI (Li)  detector 
entrance  window,  much  of  the  sample  was  masked  by  the  detector  housing. 
This  made  the  original  assumption  of  an  even  distribution  of  activity 
over  the  deposition  area  very  critical.  The  efficiencies  in  this  case 
had  relative  standard  deviations  as  high  as  13f».  It  was  felt  that  mask- 
ing off  the  sample  for  alpha  detection  to  the  same  central  7mm  diameter 
area  seen  by  the  Si (Li)  detector,  would  reduce  the  Impact  of  uneven 
deposition  on  these  calculated  ratios.  Rather  than  calculate  efficien- 
cies for  both  detectors  with  this  new  geometry,  a total  efficiency  ratio 

240 

(en»,/^n)  was  calculated  using  a 0.1  yUCi  high  purity  Pu  source  of 
the  same  geometry. 

Singles  Spectroscopy  Results 

Using  the  masked  source  for  alpha  spectroscopy,  the  calculated  L 
x-ray  to  alpha  emission  ratios  (absolute  intensities)  are  tabulated  in 
Table  V for  several  samples  of  different  Isotopic  content  and  activity. 
The  uncertainties  are  the  standard  deviations  and  incorporate  the  statis- 
tical uncertainties  in  the  total  count. 

For  the  purpose  of  demonstrating  the  relative  sensitivity  of  L 

x-ray  spectroscopy  to  gamma-ray  spectroscopy,  Table  VI  lists  the  minimum 

279  240  , 

activities  of  ^ Pu  or  Pu  necessary  for  quantitative  analysis  (for  a 

relative  uncertainty  of  1C$)  at  several  acquisition  times.  The  back- 
ground was  dependent  upon  the  activity  of  the  sample  and  this  was  in- 
corporated in  the  calculation  using  the  detector  efficiencies  for  this 
experiment. 
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Table  V 


Measured  L X-Ray  Intensities 


L x-ray/a. 

L x-ray/oc 

Sanule  # 

(measured') 

(actual  Ref:4) 

% Error 

Pu-Ol  2itO 

(800  DPM  Pu) 

! 0*0385+.  0028 

0.0378+. 0006 

+ 1.9 

L^a  s 0 . 0409+ . 0029 

0.0484+.0007 

-15.5 

L/:0.0104+.0019 

0.0120+.0003 

-13.3 

EL  j0.0898+.0045 

0.1006+. 0010 

-10.7 

Pu-05  2'iq 

(1500  DPM  Jy?u) 

L*  j 0.01 6l+. 0011 

0.01 82+. 0004 

-12.1 

Ip  :0. 021 1+, 0012 

0. 021 6+. 0004 

- 2.3 

Ly  s0.0045+.0006 

0.0053+.0001 

-15.1 

EL  :0.0417+,0017 

0.0463+.0006 

- 9.9 

Pu— 11  ?un 

(400  DPM  Pu) 

L^j 0.0438+, 0036 

0.0378+.0006 

+15.8 

Lp j0.0501+,0038 

0.0484+. 0007 

+ 3.3 

Ly  s 0.0129+. 0028 

0.0120+.0003 

+ 7.5 

EL  i0.1068+.0059 

0.1006+.0010 

+ 6.2 

Pu-13 

(1300  DPM  240Pu) 

1^  1 0,0382+, 002 6 

0.0378+. 0006 

+ 0.5 

L^sO.0521+,0031 

0.0484+. 0007 

+ 7.4 

Ly  t0.0139+. 0019 

0.0120+.0003 

+15.8 

EL  t0.1042+.0045 

0.1006+.0010 

+ 3.6 

Am-241  2, 

(650  DPM  AAm) 

I^i0.1343+.  0070 

0.132  +.003 

+ 1.5 

L^i 0.1 725+. 0074 

0.1 925+. 006 

- 9.9 

Lyi 0.0316+. 003 8 

0.0485+.002 

-3^.8 

ZL  *0.338  +.014 

0.382  +.007 

-11.5 

Table  VI 


Minimum  Sample  Activities  for  Quantification 
using  L X-Ray  Spectroscopy 


Acquisition 
Time  (Hrs) 

239Pu  (DPM) 

240Pu  (DPM) 

1 

29,409 

13,414 

4 

7,352 

3,354 

8 

3,676 

1,677 

12 

2,451 

1,118 

16 

1,838 

838 

24 

1,225 

559 

48 

612 
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Conclusions 

Coincidence  Spectroscopy.  At  the  activity  level  of  several  hundred 

239  240 

to  several  thousand  disintegrations  per  minute  of  Pu  or  Pu,  coinci- 
dence spectroscopy  appears  unfeasible.  This  is  primarily  due  to  the  low 
abundance  of  gamma  rays  from  the  decay  of  these  two  isotopes.  For  alpha 
gated  gamma-ray  spectroscopy,  the  reduction  in  background  possible  in 
coincidence  spectroscopy  is  more  than  offset  by  the  inclusion  of  the  ef- 
ficiency of  the  alpha  particle  detector.  This  reduces  the  total  number 
of  coincident  counts  and  increases  the  minimum  activity  for  quantitative 

analysis.  3ased  upon  the  background  of  the  "singles''  gamma -ray  spectrum, 

239 

the  minimum  activity  of  Pu  necessary  to  quantify  the  52  keV  gamma  ray 

in  one  hour  of  acquisition  is  only  174,000  DPM  (less  than  half  that  for 
coincidence  spectroscopy).  It  is  doubtful  that  increases  in  efficiency 
by  larger  geometries  and  higher  intrinsic  efficiencies  could  lower  these 
minimum  activities  for  coincidence  spectroscopy  by  more  than  an  order  of 
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magnitude.  These  improvements,  if  possible,  would  still  do  little  to 

improve  the  utility  of  the  system  for  analysis  of  low-activity  samples. 

239  240 

Singles  Spectroscopy.  The  principle  of  determining  Pu  to  Pu 
ratios  by  the  different  L x-ray  to  alpha  yields  is,  in  theory,  possible. 
The  calculations  are  very  sensitive  to  errors  in  the  calculated  L x-ray 
to  alpha  ratio.  It  should  be  possible  to  reduce  these  errors  to  some- 
where below  % with  better  sample-detector  geometry  for  x-ray  detection. 

With  a standard  deviation  error  of  2?o  in  the  total  L x-ray  to  alpha  ratio 

239  <240 

(Ref  4j580)  the  standard  deviation  of  the  calculated  ^ Pu/  Pu  ratio 
should  be  under  for  a 50?o  activity  contribution  of  each  isotope.  In 
conclusion,  the  sensitivity  of  the  method  is  not  as  accurate  as  mass 
spectrometry  (uncertainties  are  under  for  sample  activities  as  low  as 
100  DPM) j however,  it  does  provide  an  alternative  method  for  computing 
the  mass  ratio  of  the  two  plutonium  isotopes. 

Recommendations 

Results  of  the  "singles"  alpha  and  L x-ray  spectroscopy  approach 
have  identified  one  primary  improvement  that  is  necessary  in  reducing  the 
statistical  and  systematic  errors  in  the  measured  L x-ray/alpha  emission 
ratio.  The  improvement  would  be  that  of  increasing  the  sample-detector 
geometry.  A larger  detector  would  increase  the  geometry  factor  and  would, 
at  the  same  time,  reduce  the  sensitivity  to  irregularities  in  sample  elec- 
trodeposition. The  acquisition  times  would  be  lower  and  errors  resulting 
from  the  counting  statistics  would  be  reduced 

Two  semiconductor  detector  types  are  routinely  used  in  x-ray  spectro- 
scopy t silicon,  Si(Li),  and  germanium,  Ge(Li)  or  intrinsic  Ge.  Each  have 
their  merits  and  demerits  and  each  can  be  obtained  in  large  surface  areas 
(100-1000  nun2). 


41 


Si(Li)  detectors  exhibit  Intrinsic  efficiency  curves  that  are 
essentially  level  at  energies  from  1 keV  - 20  keV,  Above  20  keV,  the 
efficiency  rapidly  decreases  with  increasing  energy.  The  energy  at 
which  the  efficiency  starts  to  decrease  and  its  rate  of  decrease  is 
dependent  upon  detector  thickness.  The  silicon  K x-ray  escape  probabil- 
ity is  very  small,  at  less  than  1%  for  photons  of  1-40  keV  (Ref  16:2.12), 
A major  disadvantage  in  the  Si(Li)  detector  is  that  as  the  surface  area 

of  the  detector  increases,  the  resolution  becomes  much  poorer  than  that 

2 

of  a comparable  size  germanium  detector.  The  resolution  for  the  30  mm 

detector  used  in  this  experiment  was  205  eV  FcfKM  at  6.4  keV.  For  a new 

2 

Si(Li)  detector  of  500  mm  , a resolution  of  only  500  eV  FWHN  at  6,4  keV 
can  be  expected. 

Germanium  detectors,  both  intrinsic  and  lithium-drifted,  exhibit 
efficiency  curves  that  fall  off  slowly  on  both  ends  of  the  energy  spec- 
trum, Below  20  keV,  the  efficiency  falls  with  decreasing  energy  to  a 
discontinuity  at  11,1  keVj  this  is  due  to  the  K absorption  edge  of  Ge, 
With  increasing  energy,  the  efficiency  remains  relatively  level  before 
falling  off  slowly  at  60  keV.  The  rate  of  this  decrease  in  efficiency 
is,  as  with  the  Si(Li)  detector,  dependent  upon  the  thickness  of  the  de- 
tector, The  germanium  K x-ray  escape,  however,  cannot  be  neglected  as 
in  Si(Li)  detectors,  and  is  as  high  as  12f5  for  a 15  keV  incident  photon 
(Ref  16:2,12),  Germanium  detectors  have  better  resolution  than  Si(Li) 

as  the  detector  area  increases.  3emis  et  al,  (Ref  4:580)  report  a res- 

2 

olution  of  250  eV  FWHM  at  18  keV  for  a 1000  mm  intrinsic  germanium  de- 
tector. (see  Fig.  10  for  typical  Si(Li)  and  Ge  efficiency  curves). 

From  the  preceding  discussion,  it  is  obvious  that  a trade-off  in 
resolution  occurs  as  the  surface  area  is  increased.  Nearly  all 
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semiconductor-detector  L x-ray  spectra  measured  at  an  energy  resolution 
of  better  than  1 keV  FVHM  can  be  broken  down  in  four  broad  groups t the 
and  L;  groups,  even  though  all  L-series  lines  may  not  be 
readily  Identifiable.  To  identify  the  contribution  of  Np  L x rays  to  a 
spectrum  of  U L x rays,  however,  requires  a resolution  better  than  500 
eV  FWHH  at  approximately  18  keV. 

Based  upon  these  considerations,  it  would  appear  that  cost  would  be 

the  limiting  factor  on  choice  of  detector  size.  The  intrinsic  germanium 
2 

detector  of  1000  mm  described  by  Semis  would  provide  for  the  best  geome- 

2 

try,  and  its  resolution  would  be  comparable  to  the  30  mm  Si(Li)  detec- 

2 

tors.  A high  quality  Si(Li)  detector  of  200  mm  surface  area  should  how- 
ever, provide  sufficient  resolution  at  a much  lower  cost,  and  yet  in- 
crease the  geometry  greatly.  Acquisition  times  or  minimum  activities 

for  quantification  could  be  reduced  by  one  order  of  magnitude  by  an  in- 

2 

crease  in  detector  area  to  100  mm  . Table  VII  reflects  the  minimum  activ- 
ity levels  for  quantification  for  a four-hour  acquisition  time  that 
might  be  expected  for  Si(Li)  detectors  of  various  areas.  For  the  purpose 
of  table  generation,  the  background  was  assumed  to  be  similar  to  that  of 
the  detector  used  in  this  project  yet  was  assumed  to  increase  in  propor- 
tion to  surface  area  increases.  The  intrinsic  photopeak  efficiency  was 
assumed  to  be  unchanged  from  this  experiment.  Geometry  factors  for  the 
detectors  were  calculated  on  the  assumption  that  the  entire  sample  could 
be  seen  by  the  detector  and  that  the  sample-detector  distance  was  fixed 
at  5 nun. 
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Table  VII 


Minimum  Sample  Activities  for  Quantification  using  L X-Ray 
Spectroscopy  for  Various  Detector  Areas  in  a 
Four  Hour  Acquisition  Period 


2 

Detector  Area  (mm  ) 

239Pu  (DPK) 

240Pu  (DPM) 

80 

622 

285 

100 

612 

280 

200 

495 

225 

300 

380 

174 

500 

295 

135 

800 

253 

115 

1000 

240 

109 

It  is  recommended  that  work  be  continued  in  analyzing  the  L x-ray 
spectra  of  mixed  plutonium  isotopes  and  in  minimizing  the  errors  in  meas- 
ured L x-ray  Intensities.  A concentrated  effort  should  be  made  toward 
establishing  a rigid,  reproducible  source-detector  geometry,  with  dis- 
tance from  source  to  detector  determined  to  a high  degree  of  accuracy. 

A chamber  similar  to  that  used  in  this  project,  but  perhaps  with  a larger, 
thinner  beryllium  window,  to  accommodate  the  larger  geometry,  could  pro- 
vide such  a sample  mount.  A cryogenically-cooled  x-ray  detector  in  a 
second  vacuum  chamber  with  similar  beryllium  window  would  allow  simul- 
taneous alpha  and  x-ray  spectra  acquisition  (on  an  appropriate  MCA)  with- 
out the  inconvenience  of  warming  the  detector  and  reevacuating  the  cham- 
ber when  a sample  was  to  be  changed.  Emphasis  needs  to  be  continued  on 
this  simultaneous  acquisition  if  the  procedure  is  to  be  reduced  to  that 
of  a single  measurement. 
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A computer  program  similar  to  that  of  SFECTRA  should  be  developed 
to  analyze  the  resulting  complex  x-ray  spectra  resulting  from  the  pre- 
sence of  other  heavy  elements. 

Research  for  this  project  has  exposed  the  fact  that  little  is  pre- 
sently known  about  the  intensities  of  the  L- series  lines  in  each  sub- 

236 

group  of  N’p  characteristic  x rays  following  the  decay  of  Pu.  Further 
work  on  this  project  might  well  provide  information  on  these  intensities. 
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Appendix  At 

Heavy  Elegant  Radiation  Yields 


Table  VIII 


Major  Alpha  and  Photon  Energies  and  Intensities 


for  Heavy  Elements 


236Pu  1 1.-2 . 85  y 
2 

* 

Radiation 

Enerscy 

Absolute  Yield 

(Ref  1) 

5.769  MeV 

0.69 

00 

5.722 

0.31 

*157 

5.616 

0.0018 

47.  keV 

3.1X10"4  (a) 

*2 

110. 

1.2x1 0'4 

165. 

6.6x10  ^ 

238Pu  ti-  86  v 

2 

*0 

5.499  MeV 

0.72 

*43 

5.456 

0.28 

*•143 

5.358 

9. Oxl O'4 

I*  (U)  (Ref  4) 

11.618  keV 

0.0026 

1*00 

13.615 

0.0415 

00 

17.220 

0.0561 

L i (U) 

20.167 

0.0136 

h 

43.49 

3.9xl0'4  (a) 

99.87 

7.4x1  O'5 

^3  • 

152.77 

l.OlxlO"5 
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Table  VIII  (coni.) 


242  5 

Pu  ^>“3 . 79x1 0"’  y 

2 


Radiation 

Energy 

Absolute  Yield 

&0  (Ref  1) 

4.903  MeV 

0.76 

°^45 

4.863 

0.24 

Ljt  (U)  (Ref  4) 

11.618  keV 

.,0.0021 

w (u) 

13.615 

0.0310 

L/s  (U) 

17.220 

0.0415 

W (U) 

20.167 

0.0108 

^1 

44.92 

3.3X10"4  (a) 

103.50 

7.2xl0-5 

*3 

158.80 

4. 2x1 O'6 

24\m  t_i=7.95xl0^  y 

0(?5  (Ref  1) 

5.276  MeV 

0.879 

**118 

5.234 

0,106 

°^1?2 

5.181 

0.011 

°^0 

5.350 

0.0016 

°(-31 

5.321 

0.0012 

4 

74.67  keV 

0,660  (a) 

43.53 

0.055 

* 

86.7 

0.0037 

i 

117.8 

0.0050 

i 

142.0 

0.0013 

50 


Table  VIII  (cont.) 


24*Am  ti=458  y 

Radiation 

Enerscy 

Absolute  Yield 

(HSf  ° 

5.486  MeV 

0.86 

°4o3 

5.443 

0.127 

*159 

5.389 

0.0139 

*0 

5.545 

0.0025 

*33 

5.513 

0.0012 

L,  (Sp)  (Ref  6> 

11.898  keV 

0,0086 

L*(Np) 

13.944 

0.132 

V(Np) 

17.750 

0.1925 

L^(Np) 

20.785 

0.0485 

h 

59.54 

0.359 

A 

99.00 

2, 02x1 0-4 

GO 

Oo 

103.00 

1.95x1 O'4 

^9 

123.01 

1.02xi0"5 

Ao 

125.30 

3. 95x1 O'5 

246Cm  t;-5.6xl03  y 

Oi  <Ref  O 
0 

5.386  KeV 

0.81 

*45 

5.342 

0.19 

^(Pu)  <Ref  4> 

12.124  keV 

0.0021 

L*(Pu) 

14.279 

0.0333 

^(Pu) 

18.294 

0.0371 

l;(pu) 

21.417 

0.0086 

j/ 

44.55 

unknown  (a) 
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Table  VIII  (cont.) 


2kh  , 

Cm  t i=l 7. 6 y 
2 

Radiation 

Energy 

Absolute  Yield 

06  (Ref  1) 
0 

5.806  KeV 

0.767 

*43 

5.763 

0.233 

^42 

5.666 

2.3xl0"4 

L,(Pu)  <Ref  4> 

• 12.124 

0,0025 

Lr(Pu) 

14.279 

0.0386 

L^Pu) 

18.294 

0.0430 

L^(Pu) 

21.41? 

0.0103 

*1 

42.82 

2 , 0x1 0”4  (a) 

>2 

98.86 

1.3xlO"5 

*3 

152.63 

1.5X10"5 

(a)  Argonne  National  Laboratory  "Gamma-Rays 
library  listing 

by  Nuclides"  computer 
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Appendix  B» 
Figures 
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PREAMP  LINEAR  AMP  DELAY  AMP  BIASED  AMP 


Figure  6.  Block  Diagram  of  X-Ray  Spectroscopy  Circuit 


J 


PHOTON  ENERGY  (keV) 


Figure 
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* MRC  ^Am  standard 
241 

H NEN  AAm  standard 

8.  Intrinsic  Photopeak  Efficiency  Curve  for  Princeton 
Gamma -Tech  Si (Li)  Detector 


6i 


L 2 3 4 5 6 7 8 > 10  20  30  40  50  (.0  100 

PHOTON  ENERGY  (keV) 

* Refj  ORTEG  Technical  Data  sheet  Series  ' " ‘0  .1(1  l)  X Hay  De tectors , Dec  19?1 

**  Hef 1 ORTEG  Technical  Data  sheet  Serle..  • .e(Ll)  Low-Energy  Photon  Detectors,  Oct  19?0 


Figure  11,  Typical  Uranium  L X-Ray  Spectrum  from 
Low-Activity  ^39pu  Sample  on  Princeton 
Gamma-Tech  Si(Li)  Detector 


5.15  5A9  5.76 


ALPHA  ENERGY  (MeV) 

Figure  12,  Typical^ 240pu  Qpha  Spectrum  with  ^^Am 
and  impurity  Traces  taken  with 

450  mnr  ORTEC  Ruggedized  Silicon  Surface 
Barrier  Detector 
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Figure  14.  Typical  Complex  U/Np  L X-ftay  Spectrum 
from  combined  ^39pu  and  24l^m  sanp]_e 
taken  with  Princeton  Gamma-Tech  Si(Li) 
Detector 
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45  59.5  104 

PHOTON  ENERGY  (keV) 


Top:  Ungated  Singles  Gamma-Rav  Spectrum 

of  combined  2“*0pu  and  24lAm  sample 
Bottoms  5.15  MeY  alpha  Gated  Gamma-Ray 

Coincidence  spectrum  of  combined 
240pu  and  ^4l^m  sample 

Figure  15.  Singles  Gamma  vs.  Coincidence  Gapma-Ray 
Spectrum  for  combined  240pu  and  ^l^m 
Sample  taken  with  ORTSC  450  mm2  Silicon 
Surface-3arrier  Detector  for  Alpha  Par- 
ticles and  Nuclear  Diode  Ge(Ll)  Detector 
for  Gamma  Rays. 
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APPENDIX  C 


GEOMETRY  FACTOR  FOR  AN  EXTENDED  DISK 
SOURCE  ON  AXIS  WITH  A DISK  DETECTOR 

Many  of  the  source-detector  geometries  with  which  this  project 
dealt,  could  not  be  accurately  approximated  by  point  source  geometries. 

A simple  program  that  would  calculate  this  geometry  factor  given  the 
radii  of  detector  and  source  and  the  distance  between  them,  proved  to  be 
of  value  in  determining  detector  efficiencies. 

The  approach  taken  was  to  simplify  the  expression  presented  by 
Heath  (Ref  15:21)  for  the  energy-dependent  efficiency  of  a disk  source 


The  energy  dependence  of  the  equation  can  be  removed  by  setting  the 
square  bracketed  terms  under  the  integrals  equal  to  one.  These  are  the 
probabilities  (as  a function  of  distance)  that  the  radiation  will  be  ab- 
sorbed. Then,  by  taking  the  limit  at  t -*■  0 (for  an  infinitely  thin 
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detector),  equation  (12)  simplifies  to  the  geometry  factor  (GF)  effi- 
ciency 


TT 

■2 


GF  = 


itr* 


n 

r r 

xdx 

d^> 

J 

_i  ( - xsiw <t>  +■  /xJ5UAlqS  — ( x1- ) 


stwe-d©- 


(13) 


Integrating  over  0 the  triple  integral  is  reduced  to  the  double  integral: 


Re 

r 


GF  = 


1TR-* 


J 


2 

r 


xdx 


I — Cos 


(1^) 


simplifying  this  integral  provides : 


r 


i 


<S|F  = 


I 

2. 


TTR^ 


Xdx 


Cos 


w/ 

0 


J 

_TT 

2 


to N ^ (X-Rj)  ||d0  (15) 

k """ 


This  double  integral  was  then  solved  by  calling  on  AFIT  Subroutine  SIMPD 
which  solves  double  integrals  by  Simpson's  Rule  iteration.  The  ource 
program  follows  in  this  Appendix. 
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PROGRAM  SOL  I DA  ( INF'UT-=/30 . OUTPUT ) 

COMMON  /l/  R D > H 

THIS  PROGRAM  COMPUTES  THE  GEOMETRY  FACTOR  FOR  A 
CIRCULAR  DISK  SOURCE  A GIVEN  DISTANCE  AWAY  AND 
ON  AXIS  WITH  A CIRCULAR  DISK  DETECTOR. 

where:  r s = s n u r c f.  r a d i u s 

RD- DETECT OR  RADIUS 

H=DISTANCE  FROM  SOURCE  TO  DETECTOR 
PI =3. 1415926535897932384a 
5 READ#  » RS » RO  * H 

IF  (EOF (5L INPUT)  ) 99.10.99 
10  PRINT#. "SOURCE  RADIUS  IN  MH='»RS 
PRINT*. 'DETECTOR- RADIUS  IN  MM-'.RD 

PRINT*. 'DISTANCE  FROM  SOURCE  TO  DETECTOR  IN  MM-'.H 
A 1 ~ - P I / 2 . 

Bl=PI/2 . 

A2  = 0 
B2  = RS 

CALL  SIMPDCAJ.  .EM  .A2.B2.  .001 . .001 . 100 . 100 . 2 . 2 . ANSQ . 
1 ALARM  » JTOL ) 

PRINT*. ' ALARM- ’ .ALARM.  " JTOL~* .JTOL 
SOA-1 . /2 . -ANSQ/ ( P I*RS**2 ) 

PRINT*. "THE  GEOMETRY  FACTOR  TSi'.SOA 
20  PR I N T * » ' * * * * * * ********  * * * * * * * * * * * * * * * * * * * * * * * * * * " 
GO  TO  5 
99  STOP 
END 


SUBROUTINE  FNEVAL  (PHI.FPHIX) 

COMMON  X 

COMMON  /!/  RD.H 

FPHIX  = X*COS(ATAN( <-X#SIN<PHI ) 

8+ SORT (AB3(X**2*(SIN(PHI) ) **2- ( X**2-RD**2 ) ) ) )/H) ) 
RETURN 
END 
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APPENDIX  D 


SPECTRA  - A PROGRAM  TO  PRODUCE  A THEORETICAL  L X-RAY  SPECTRUM 
FROM  THE  SUM  OF  GAUSSIAN  FSAXS  AND  LEAST-SQUARES  FIT  BACKGROUND 

Program  SPECTRA  was  written  to  provide  a method  of  correcting  a 
complex  L x-ray  spectrum  for  the  presence  of  other  nuclides.  Given  a 
magnitude  and  energy  of  a reference  peak,  the  program  takes  the  energies 
and  relative  intensities  of  a series  of  L-shell  lines;  fits  Gaussian 
peaks  to  these  x rays  and  sums  them  to  provide  a theoretical  spectrum. 
Given  the  background  count  at  respective  channel  numbers,  SPECTRA  calls 
on  a third  degree  polynomial  least-squares  fit  to  the  background  points 
(AFIT  Subroutine  PLSQ)jadds  the  background  to  the  theoretical  spectrum, 
then  subtracts  that  total  spectrum  from  an  experimental  spectrum  read 
into  the  program  on  data  cards.  The  resulting  spectrum,  then,  consists 
of  the  peaks  of  interest.  Given  areas  of  interest,  SPECTRA  then  sums  the 
counts  under  the  resulting  peaks. 

FUNCTION  PGAUSS  (Ref  12s45)  provides  the  Gaussian  curve  as  a func- 
tion of  energy  and  standard  deviation.  The  standard  deviation  also  de- 
pends upon  the  energy  of  the  peak  and  is  provided  to  PGAUSS  through 
FUNCTION  SIG.  The  standard  deviation  is  proportional  to  the  square  root 
of  the  energy.  The  constant  of  proportionality  in  SIG  is  supplied  by  the 
user  and  is  dependent  upon  total  detector  system  resolution. 

Fig,  16  depicts  a theoretical  Np  L x-ray  spectrum  and  background 

241 

arising  from  the  decay  of  Am.  The  relative  intensities  and  energies 
of  the  16  L-shell  lines  used  to  construct  this  spectrum  were  taken  from 
Hyde  et  al.  (Ref  11 *853) • Fig.  17  shows  this  same  spectrum  fit  to  the 
Np  L x-ray  peaks  in  a complex  U/ttp  L x-ray  spectrum.  Fig.  18  provides  the 
net  U L x-ray  spectrum  after  subtracting  the  background  and  Np  spectrum. 
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P R n G R A M S F' E: C T R a ( I N P U T =- / 8 0 , DUT P LIT  , P L Q T ) 

THIS  PROGRAM  COMPUTES  AND  PLOTS  A GAUSSIAN  DT SIR [BUTTON 
FOR  AN  ENERGY  PEAK  GIVEN  THE  MEAN  ENERGY.  AND  RELATIVE 
INTENSITY  OF  EACH  PEAK  IN  THE  SUBGROUP  AND  THE  MAGNITUDE 
OF  THE  HIGHEST  PEAK  IN  THE  SUBGROUP.  THE  PROGRAM  THEN 
SUMS  THE  INDIVIDUAL  GAUSS I AN 3 AND  PLOTS  THE  TOTAL  SPECTRUM 
D I M E N S 1. 0 N G A ( 1 0 3 0 ) , X C ( J.  0 3 0 ) - S P E C ( 1 0 30).  B K G S ( 1 030) 

D I M E N S 1 0 N X B ( 9 ) . Y B ( 9 ) » C B ( 4 ) 

D I M E N S I ON  S P E CIN(  1. 0 3 0 ) . S S P E C ( J.  0 3 0 ) 

1 F 0 R M A T ( 1 H . 1.0  ( 1 P E 1 3 . -1 ) / 1. 0 ( :L  P E13 . 4 ) ) 

2 FORMA  T(.1H  - F6 . 2 . 5X , 1PE1 3 . 4 ) 

8 FORMAT  < 1F6.0,9FQ.0> 

9 FORMAT  (III  , t OF  13.0  /10(F13.0)) 

SLOPE-. 01 36479758 
YINCPT=10. 24680  13 

GA  ( 1025  ) -'I  . $ SPEC  ( 1020)  = 1 . $ BKGS  < 1025 ) = 1 . 

GA(1026)  = 1.  * SPEC ( 1026 ) = 1 . $ BKGS ( 1026  ) =1  . 

READ  8. (SPEC IN ( I ) ,1=1.1024) 

DO  91  1=1,1023 
91  SPEC  IN  ( I ) =SF’ECIN  ( I + 1 ) 

READ* , N 
DO  3 1 = 1,  N 

3 READ* , XB ( I ) , YB ( I ) 

CALL  PLSQ ( XB , YB , N , 3 , CB , 0 , EMAX , ERMS , EMEU ) 

DO  4 1=1,1024 
X=FLQAT ( I ) 

XC  ( I ) =FL.OAT  ( I ) 

4 BKGS ( I > =CB ( 1 )*X**3+CB(2)*X**2+CB(3)*X  + CB<4> 

5 READ* , NPKS  » F’KMAG 

IF ( EOF (5L INPUT) ) 199, 7, 199 
7 DO  90  L=1 , NPKS 
READ* , EAVG  » R I 

PRINT*. ***********  PEAK  NUMBER  *,N+1,*  DATA  ********** 
PRINT*.*  PEAK  ENERGY** , EAVG, * KEV* 

PRINT*. * PEAK  RELATIVE  INTENSITY3* ,RI 

PRINT*,*  PEAK  MAGNITUDE®* ,PKMAG*RI , * COUNTS’ 

S I G=S I GMA ( EAVG ) 

PRINT*. * S I GMA= " , SI G 

PRINT*.  * FWHM = * , 2 . 35*SIG  , *KEV 

SUM  = 0 

DO  10  1=1,100 
X=EAVG-3.*SIG+.06*I*SIG 
RJ=(X-YINCPT) /SLOPE 
J=R  J 

• IF ( < R J-J ) . GE . . 5 ) J*R J+ 1 

10  GA ( J ) =PKMAG*R I *SIG*PGAUSS ( X , EAVG , S IG ) / . 3989422904 
DO  20  1=1,1023 

20  SUM--SUM+  ( GA  ( I ) +GA  ( 1 + 1 ) >/2. 


PRINT* • * TOTAL  COUNTS  UNDER  GAUSSIAN* ' » SUM 
IF  (K.GF.l)  GO  TO  AO 
CAl  L PLOT ( 4 » rt.r~3) 

CALL  SCAL  F ( X C 7 6 . 0 * J.  02 4 » 1 ) 

CALL  AX I S < 0 . . 0 . f 1 4 H CHANNEL  NUMBER »-14.6.0»0» XC ( 1 025 > 
XrXC ( l 0 2 6 > ) 

CALL  LG  AX  I S < 0 . . 0 . 7 1. 6HC0UNTS  / CHANNEL  7 1. 6 ? 5 . . 90 . r 
SGA ( 1025) 7 GA ( 1026) ) 

*0  DO  25  [ -17  1024 

SPEC ( I ) -GA ( I ) + SPEC  ( I ) 

I F ( £< K G 3 ( C ) . LE . 1 . ) B K G S ( I ) - J.  • 

25  IF  ( GA  < I ) . LE . 1. , ) GA  ( I ) -1  . 

CALL  LGL I NE  ( XC  7 GA  7 .1 02 4 r 0 7 1 7 1 ) 

DO  50  J-1.7  1024 
30  G A ( J ) = 0 
K-K+l 

80  CONTINUE 

99  PRINT* 7 'THE  TOTAL  SPECTRUM' 

DO  40  J=:  1 7 1 024 

SPEC ( J ) =SPEC ( J ) +BKGS ( J ) 

40  SSPEC  < J ) --SPEC I N ( J ) -SPEC  ( J ) 

PRINT  If  SPEC  ( J ) 7 J—  1 71024) 

READ *7  NSEC 

DO  100  L ~ 1 7 M S E C 

SECS*0 

RSECS=0 

READ*7N1 7N2 

DO  101  K=N1 7 N2 

RSECS-RSECS  + SSPEC ( K ) 

101  SECS=SECS+SPEC(K) 

PR  I NT*  7 “ SUM  COUNTS  CH.VN^'  TO  CH.'7N27*=  *7SECS 
PR INT*  7 ' THE  REMAINING  SUM  COUNTS  AFTER  STRIPPING  * 

100  PRINT*. •CHANNEL" 7N1 7 " TO  CHANNEL " 7 N2 7 1 = * 7RSECS 
CALL  LGL I NE ( XC  7 BKGS 7 1 024  7 - 1 00  7 4 7 1 ) 

CALL  LGLINE(XC7SPEC  7 102  4 7 -07  3 7 1 ) 

CALL  PLOT (9. , 0. 7-3) 

87  DO  88  1=1 7 1024 

IF  ( SSPEC ( I ) . LE . 1 . ) SSPEC ( I ) = 1 . 

88  IF < SPEC IN ( I ) . LE . 1 . > SPECIN(I)  = 1. 

PRINT  9 7 ( SPEC  I N ( I ) 7 1 = 1007200) 

SPECIN( 1025)  = 1 . 

SSPEC < 1025)  = 1 . 

SPEC IN ( 1026 ) = 1 . 

SSPEC ( 1026 ) =1 . 

CALL  AXIS(0. 7 0 . fl 4HCHANNEL  NUMBER  7 -1 4 7 6 . 0 7 0 7 XC ( 1025 ) 7 
SXC< 1026) ) 

CALL  LGAX IS  < 0 . 7 0 . 7 16HC0UNTS  / CHANNEL  7 1 6 7 5 . 7 90 . 7 
SGA ( 1025) 7 GA( 1026) ) 
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CALL  LGl.  1 NE  < XC.  , SPEC  , 1 024  ,-0,3,1) 

CALL  LGLI  NE ( XC  r SPEC  IN  * 1 024  > -2  > 3 * 1 ) 

CALL  PLOT ( 9 . ,0. *-3> 

CALL  AX  I S ( 0 . • 0 . * 1 AHCI-I  -iNNEL  NUMBER  » - 1 4 » 6 . 0 , 0 • XC  ( 1 025  > * 
aXC (1026) ) 

CALL  LG AX IS ( 0 . , 0 . > 1 6MC0UNTS  / CHANNEL » 1 6 ? 5 . * 90 . » 

*GA(  1025)  , GA  ( .1.026:  ) 

CALL  LGLTNE < XC  , SSPE C , 1024 >-0,3, 1 ) 

CALL  PLOT (9. ,0. * -3) 

CALL  PLOTE 
19?  STOP 
END 


FUNCTION  P G A U S S < X r E A V G r S I G ) 

DOUBLE  PRECISION  Z 
Z=<  X-EAUG ) /STG 

PGA US S- . 39S?422804/SIG*DEXP ( - ( Z**2 >/2.  > 

RETURN 

END 


FUNCTION  SIGMA  < EAUG ) 

F~EXP ( . 0607 1KEAUG-3. 081 16) 

SIGMA=3 . 6*S0RT < . 0029*EAUG*F ) /2 . 35 

RETURN 

END 
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COUNTS  / CHANNEL 
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Figure  1?.  Thforetical  Fit  of  Np  L X-Ray  Spectrum  to  Experimental 
Complex  Np/U  Spectra 
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COUNTS  / CHRNNEL 


Figure  18.  Uranium  L X-Ray  Spectrum  after  Stripping  away  Neptunium 
L X-Ray  Spectrum  from  Complex  Spectra 
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